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ABSTRACT 


This  is  the  final  report  covering  research  directed 
toward  the  study  of  the  seismicity  of  the  Southeastern 
United  States.-^  Travel-times  determined  from  local  earth¬ 
quake  and  refraction  data  are  pr&sented  which  indicate 
a  cruotal  structure  of  h^  ■  33.0  km  (a  »  5.88  km/sec), 
h^  *  10.8  km  (a  -  6.58  km/sec),  and  an  upper  mantle  velocity 
of  8.10  km/sec.  Fundamental  and  first  higher  order  Rayleigh 
group-velocity  data  determined  oy  digital  bandpass  filter¬ 
ing  are  presented  for  the  Southern  Appalachian  region.  The 
Dunkin  modification  of  the  Thomson-Haskell  matrix  method 
is  used  to  compute  theoretical  Rayleigh  dispersion  curves 
for  comparison  with  the  observed  curves.  A  slight  velocity 
reversal  in  the  upper  crust  centered  at  about  15  km,  a 
general  increase  of  crustal  velocities  and  densities  with 
depth  below  this  zone,  and  an  upper  mantle  low  velocity 
zone  beginning  at  a  depth  of  70  km  are  indicated  beneath 
the  Southern  Appalachians.  The  Appalachian  foreland  has 
crustal  structure  similar  to  the  Gutenberg-Birch  II  conti¬ 
nental  model  with  a  total  thickness  of  40  km. 

A  sin  x/x  analysis  of  the  Bouguer  gravity  data  yields 
a  total  crustal  thickness  of  about  50  km  beneath  the  Southern 
Appalachians . 

P- residuals  computed  at  Chapel  Hill,  North  Carolina  and 
McMinnville,  Tennessee  show  a  systematic  deviation  of  as 
much  as  +  3  sec. 
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1.0  INTRODUCTION 


This  final  report  covers  the  work  done  on  contract  AF  19(628)-3892, 
"Research  Directed  Toward  the  Study  of  the  Seismicity  of  the  Southeastern 
United  States".  The  report  is  mainly  concerned  with  the  work  done  over  the 
past  year;  the  determination  of  crustal  and  upper  mantle  structure  in  the 
Southeastern  United  States.  The  First  Annual  Technical  Report  [Minear,  1965] 
covers  the  development,  installation,  and  calibration  of  the  short-period 
displacement  seismograph  at  the  University  of  North  Carolina  and  the  RTI 
field  refraction  system.  The  Second  Annual  Technical  Report  [Minear,  1966] 
includes  the  location  of  local  epicenters,  the  results  of  the  field  refraction 
studies,  the  results  of  the  computation  of  P-residuals,  the  calculations  of 
magnitude,  focal  depth,  ar.d  energy  release  for  several  local  earthquakes, 
and  preliminary  crustal  structure  estimates  from  gravity  data. 

Research  accomplishments  during  the  performance  of  the  contract  are 
briefly  summarized  below. 

1)  During  the  first  year  of  work  a  short-period  displacement  seismo¬ 
graph  system  was  designed,  constructed,  and  placed  on  routine  operation  at 
the  University  of  North  Carolina  seismograph  vault  at  Chapel  Hill,  North 
Carolina.  Although  the  system  performed  well  [Minear,  1965],  the  background 
noise  level  at  the  UNC  station  was  too  high  to  permit  recording  of  local 
earthquakes  located  mainly  in  the  Southern  Appalachians.  At  present, 

a  remote  vault  is  currently  under  construction  by  the  University  to  provide 
an  up-to-date  seismic  facility. 

2)  Refraction  work  was  carried  out  using  local  quarry  blasts  as 


energy  sources. 


3)  Local  travel-time  curves  were  developed  using  several  of  the  major 
local  earthquakes  which  were  well  recorded  by  portable  and  permanent  stations 
in  the  region. 

4)  P-residuals  were  computed  for  several  hundred  epicenters  recorded 
at  the  Cumberland  Plateau  Seismoloa^cal  Observatory  and  Chapel  Hill,  North 
Carolina.  A  systematic  deviation  of  the  residuals  similar  to  that  noted 
by  other  investigators  was  found.  This  deviation  cannot  be  explained  by 
crustal  velocity  variations  and  must  indicate  a  real  error  in  the  Jeffreys- 
Bullen  travel- times. 

5)  Estimation  of  focal  depth,  magnitude,  and  energy  release  from 
previous  intensity  studies  of  four  Southeastern  earthquakes  were  made. 

6)  Total  thickness  of  the  Southern  Appalachian  crust  was  determined 
using  Bouguer  gravity  anomalies  and  the  sin  x/x  method  of  computing  the 
mass  anomaly  producing  a  given  gravity  anomaly. 

7)  Crustal  and  upper  mantle  structure  was  determined  using  fundamental 
and  first  higher  Rayleigh  mode  group  velocity  dispersion. 

8)  Computer  programs  were  written  for  bandpass  filtering,  computation 
of  P-residuals,  least  squares  epicenter  location,  computation  of  theoretical 
travel-times  from  a  given  velocity  structure,  computation  of  theoretical 
Rayleigh  dispersion  curves  and  modal  shape,  computation  of  the  variation 

of  phase  velocity  with  layer  parameter  variations,  and  the  computation  of 
the  mass  anomaly  from  a  given  gravity  anomaly  profile. 

At  the  start  of  the  project,  it  was  anticipated  to  do  considerable  work 
on  phase  and  amplitude  spectra  of  bot:'  seismic  signals  and  background  noise. 
Also,  it  was  hoped  that  more  work  could  have  been  done  on  general  seismicity, 
distribution  of  epicenters,  tocal  depth,  and  energy  release.  Failure  co 
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much  work  had  been  done  on  back- 


acquire  a  digital  system  and  the  fact  that 
ground  noise  did  not  make  the  study  of  spectra  appear  worthwhile.  The 
general  seismicity  study  was  frustrated  by  the  poor  recording  station  dis¬ 
tribution  in  the  region.  Therefore,  crustal  and  upper  mantle  structural 

dies  utilizing  re  raction.  gravity,  and  surface  wave  dispersion  were 
concentrated  on. 

This  report  specifically  covers  the  local  travel-time  curves  for  the 
Southern  Appalachian  region  (Sec.  2).  the  determination  of  crustal  thickness 
from  gravity  data  (Sec.  3) ,  the  computation  of  theoretical  dispersion  curves 
(Sec.  4),  the  determination  of  Rayleigh  group  velocities  (Sec.  5),  and  the 
crustal  and  upper  mantle  structure  in  the  Southeastern  United  States  (Sec.  6) 
Numerical  computational  methods,  tables,  charts,  and  computer  program 
listings  are  presented  in  the  Appendices. 
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2.0  TRAVEL-TIME  CURVES  FOR  THE  SOUTHEASTERN  UNITED  STATES 


Travel-time  curves  were  determined  by  using  data  from  three  local 
earthquakes  which  were  well  recorded  by  Worldwide  Standard  Seismograph 
Stations  and  portable  Vela  stations  operating  in  the  Southeastern  United 
States  [Minear,  1966],  Fig.  4  shews  the  location  of  the  epicenters  and 
recording  stations.  Travel-time  curves  drawn  from  the  local  earthquake 
data  are  shown  in  Fig.  1.  Refraction  data  obtained  from  quarry  blasts  and 
during  the  East  Coast  Onshore  Offshore  Seismic  Experiment  and  theoretical 
travel-times  computed  for  a  typical  linear  mountain  from  the  Herglotz- 
W.iechert  equations  are  also  plotted  in  Fig.  1. 

Travel-time  curves  corresponding  to  arrivals  from  the  first  crustal 
layer  and  from  the  crust  mantle  boundary  are  drawn  from  first  arrivals 
and  are  estimated  accurate  to  within!  *1  km/sec.  Second  arrivals  were  used 
to  define  curves  corresponding  to  two  major  crustal  layers.  No  major  third 
layer  in  the  crust  is  indicated  by  the  refraction  and  earthquake  data. 
However,  first  arrivals  from  the  refraction  profiles  and  second  arrivals 
from  250  to  550  km,  indicate  that  the  crustal  velocity  may  increase  r* ;her 
continuously  from  about  10  km  to  around  45  km.  The  local  travel-time  data 
yields  a  crustal  model  of  h^  *  33,0  km  (a  ■  5.88  km/sec),  h ^  m  10.8  km 
(a  ■  6.58  km/sec) ,  and  an  upper  mantle  velocity  of  8.10  km/sec.  As  can 
be  seen  from  Fig.  4,  the  epicenters  and  recording  stations  are  generally 
located  to  the  west  of  the  core  of  the  Appalachians.  Crustal  structure 
determined  from  the  travel-time  data  thus  corresponds  to  the  crust  beneath 
the  Appalachian  foreland. 

Velocity  structures  of  the  crust  and  upper  mantle  for  the  Appalachian 
foreland,  the  Northern  Alps  (N),  Central  Alps  (C) ,  and  Northern  Alpine  fore¬ 
land  (F)  [Knopoff,et  al,  1966],  and  a  linear  mountain  belt  are  shown  in 
Fig.  2.  Crustal  velocities  for  the  Appalachian  foreland  generally  agree 
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Fig.  2  COMPRESSIONAL  VELOCITY  MODELS  FOR 
LINEAR  MT.  BELT,  APPALACHIAN  FORELAND, 
AND  ALPINE  STRUCTURE  (KNOPOFF,  etAL, 1966} 


with  those  for  the  Northern  Alps  at  depths  greater  than  about  2  km.  The 
disagreement  for  depths  less  than  2  km  can  probably  be  accounted  for  by 
the  sedimentary  cover  present  in  the  Northern  Alps.  The  5.85  km/sec  layer 
in  the  Appalachian  foreland  ,  extending  to  a  depth  of  33  km,  is  thicker 
than  any  of  the  Alpine  structures.  However,  as  mentioned,  the  Appalachian 
foreland  velocities  may  increase  rather  continuously  from  about  10  to  40  km. 
Upper  mantle  depth  in  the  Appalachian  foreland  is  greater  than  beneath  the 
foreland  to  the  north  of  the  Alps  by  14  km  and  greater  than  beneath  the 
central  Alps  by  4  km. 

A  preliminary  summary  of  seismic  refraction  work  in  the  vicinity  of 
the  Cumberland  Plateau  Seismological  Observatory  [Bf  rcherdt  et  al,  1966 j 
indicates  a  crustal  model  of  h.,  ■  12  km  (v^  -  6.1  km/sec);  *  28  km  (v£  * 

6.7  km/sec)  and  an  upper  mantle  velocity  of  8.0  +  km/sec. 
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3.0  CRUSTAL  THICKNESS  FROM  GRAVITY  DATA 

Total  crustal  thickness  was  computed  from  Bouguer  gravity  values 
along  a  Northwest-Southeast  profile  extending  from  about  460  km  off  the 
North  Carolina  coast  (33°N,  73°W)  to  the  Kentucky-Illinois  border  (38°N, 

88°W).  The  sin  x/x  method  of  Tomoda  and  Aki  [1955]  was  used  to  compute 
the  depth  to  a  mass  anomaly  producing  the  observed  gravity  anomalies. 

Bouguer  gravity  values  were  taken  from  the  American  Geophysical  Union 
Bouguer  Gravity  Anomaly  Map  of  the  United  States.  Fig.  3  shows  the  gravity 
profile  values,  tie  total  crustal  thickness  computed  from  these  anomalies, 
and  regional  subsurface  geology.  The  subsurface  geological  information  was 
obtained  from  McGuire  and  Howell  [1963]  in  Kentucky,  Hersey,  et  al  [1959] 
for  the  North  Carolina  continental  margin,  and  from  the  geologic  map  of  North 
Carolina.  Crustal  structure  to  the  crust-mantle  boundary  at  location  H' , 
and  to  2  km  at  12-13  is  based  on  refraction  profiles  of  Hersey,  et  al  The 
subsurface  geology  in  North  Carolina  is  intended  only  to  indicate  possible 
near  surface  relations  between  geology  and  Bouguer  gravity  anomalies. 

In  the  sin  x/x  method,  crustal  thickness  is  computed  from 

d(nx)  *  d  -  d’(nx)  3.0-1 

where  d  is  an  assumed  thickness  and  d'(nx)  is  a  correction  to  this  thickness 
given  by 

d’(nx)  -  .  3.0-2 

M(nx)  is  the  convolution  of  the  observed  gravity  anomalies  Ag(qAx)  with  a 

symmetric  function  $  which  is  a  function  of  assumed  crustal  thickness  and 

n 

station  spacing.  Thus, 

q«m 

M(nx)-^v  ,L  48<q4x) 


8 


Oscillations  in  the  gravity  values  Aq(qAx)  due  to  near  surface  density 
irregularities  can  result  in  Aq(qAx)  going  positive  and  negative,  such  as 
in  the  region  over  the  slate  belt  in  Fig.  3.  Thus,  M(nx)  may  oscillate 
between  positive  and  negative  values  which  in  turn  yields  positive  and 
negative  oscillations  of  d'(nx).  The  ultimate  effect  is  that  in  the  regions 
of  local  near  surface  perturbations,  the  total  crustal  thickness,  d(nx) ,  may 
oscillate  widely  about  the  assumed  thickness,  d  ,  as  can  be  seen  from  3.0-1. 
Introduction  of  a  density  contrast  Ap  * hich  varies  with  depth  will  not 
eliminate  these  oscillations.  One  must  either  choose  the  station  spacing 
wide  enough  so  that  the  convolution  of  the  anomalies  with  the  function  ♦ 
effectively  filters  out  the  local  perturbations  or  smooth  the  total  depth 
function  d(nx) . 

From  Fig.  3,  it  can  be  seen  that  the  Carolina  Slate  Belt  is  associated 
with  a  gravity  high  which  effectively  introduces  a  positive  perturbation  on 
the  regionally  decreasing  gravity.  If  the  crustal  thickness  is  computed 
without  removing  this  perturbation,  the  thickness  oscillates  about  the 
assumed  thickness  beneath  the  belt.  As  shown  in  Fig.  3,  the  local  gravity 
high  over  the  slate  belt  can  be  largely  accounted  for  if  the  belt  is 
approximated  by  a  two-dimensional  block  with  lateral  extent  equal  to  the 

_3 

slate  belt,  an  8  km  depth  and  a  density  contrast  of  .26  gm/cm  . 

Crustal  thickness  was  computed  from  3.0-1  using  3.0-2  and  3.0-3  with 
an  assumed  crustal  thickness  of  45  km  and  a  station  spacing  of  60  km.  The 
thickness  values  were  then  smoothed  with  a  three  point  moving  average  filter 
(.25,  .50,  .25)  resulting  in  the  smoothed  crustal  thickness  curves  shown  in 
Fig.  3.  Two  curves  are  plotted,  corresponding  to  crustal-upper  mantle 

_3 

density  contrasts  of  .3  and  .6  gm/cm  .  Since  the  ocean's  crust  is  denser 
than  the  continental  crust,  the  &".3  curve  approximates  the  crustal  thickness 
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better  u  der  the  continental  margin.  The  agreement  with  the  thickness  as 
determined  by  refraction  work  [Hersey,  et.al.,  1959]  at  point  H'  is  good. 
Both  curves  indicate  a  crustal  thickness  of  at  least  iO  km  under  the 
Southern  Appalachians.  Perturbations  in  the  crustal  thickness  are  caused 
ly  the  Cincinnati  Arch  and  the  Carolina  Slate  Belt.  The  thinning  of  the 
crust  necessary  to  produce  the  Bouguer  anomaly  over  the  Cincinuati  arch 
and  the  slate  belt  is  about  9  km  for  the  ^".6  curve  and  about  4.5  km  for 
the  ^  ".3  curve.  Due  to  the  magnitude  of  the  crustal  thickness  changes, 
it  appears  that  the  sources  of  the  local  highs  over  the  Cincinnati  Arch 
and  the  Carolina  Slate  Belt  are  relatively  near  surface. 

The  crust  thus  thickens  from  about  83  km  at  the  North  Carolina  coast 
to  about  51  km  beneath  the  core  of  the  Appalachians  and  then  thins  to  about 
43  km  at  the  Kentucky-Illinois  border. 


4.0  DETERMINATION  OF  RAYLEIGH  GROUP  VELOCITY 


The  locations  of  permanent  recording  stations  in  the  Southeastern  United 
States  lie  principally  along  the  Appalachian  trend  (See  Fig.  4).  Permanent 
worldwide  standard  stations  capable  of  recording  long-period  seismic  signals 
are  located  at  Spring  Hill,  Alabama  (SHA);  Atlanta,  Georgia  (ATL) ;  McMinnville, 
Tennessee  (CFO)  ;  Blacksburg,  Virginia  (BLA) ;  and  Oxford,  Mississippi  (OXF) . 
Portable  long-period  units  have  been  operated  by  the  Geotechnical  Corporation 
under  project  Vela,  but  these  stations  are  also  located  along  the  Appalachians. 
Because  of  the  widely  space  station,  it  was  impossible  to  calculate 
phase  velocities  directly  using  triangular  arrays  of  stations.  Therefore, 
epicenters  were  selected  to  give  travel  paths  parallel  or  perpendicular  to 
the  Appalachian  trend.  It  was  hoped  that  variations  in  crustal  and  upper 
mantle  structure  between  stations  located  along  the  Appalachians  could  be 
detected  by  observing  the  variation  of  group  velocity  of  a  wave  train  travel¬ 
ing  the  scation  sequence  SHA  -ATl  -BLA  parallel  to  the  Appalachian  trend  or  by 
comparing  group  velocities  at  BLA  and  ATL  from  waves  arriving  perpendicular 
to  the  Appalachian  with  those  of  a  normal  continental  structure. 

On  the  basis  of  epicentral  location,  signal  amplitude, and  availability 
of  records,  two  epicenters  were  selected  for  study.  Table  II  gives  the 

information  pertinent  to  these  epicenters. 

Table  II 


Epicenter 

Location 

(USGS) 

Date 

Time 

(USGS) 

Magnitude 

(USGS) 

Focal 

Depth(km) 

(USGS) 

Distance (km) 
to 

station 

Jalisco,  Mex. 
17. 8N,  105. 9W 

11  Oct.,  1963 

-  - 

10:17:07.6 

5.0 

33 

BLA  -  3286.3 
OXF  -  2467.3 
SHA  -  2291. C 
ATL  -  2757.6 

jS.  Alaska 
j  62 . 7N,  152. 0W 

29  June,  1964 

07:21:32.8 

5.6 

33 

BLA  -  5488.0 
OXF  -  5276.3 
SHA  -  5110.0 
ATL  -  5633.9 
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LOCATION  MAP  OF  PERMANENT  AND  TEMPORARY  SEISMOGRAPH  STATIONS 
USED  IN  TRAVEL -TIME  AND  RAYLEIGH  WAVE  DISPERSION  STUDIES 

Fig.  4 


Records  from  the  standard  stations  were  hand  digitized  at  two  second 
intervals  using  a  plastic  grid  overlay.  This  digital  data  war.  stored  on 
magnetic  tape  for  processing.  Since  the  azimuths  of  the  epicenters  measured 
from  the  recording  stations  did  not  coincide  with  either  of  the  horizontal 
component  seismograph  orientations,  Rayleigh  wave  motion  was  contaminated 
by  Love  wave  motion.  In  order  to  separate  the  Rayleigh  and  Love  wave 
motions , radial  and  transverse  seismograms  were  generated  from  the  North- 
South  and  East-West  components  at  each  station.  The  relations  used  in  the 
transformations  were 

radial  component  ■  r  ■  OE  cos 6  +  ON  sin0  ,  and 

transverse  component  «  t  ■  OE  8in0  -  ON  cos6  . 

where 

6  *  azimuth  of  epicenter  from  station  measured  counterclockwise 

from  east, 

OE  =  east-west  component  amplitude,  positive  toward  east, 

ON  ■  north-south  component  amplitude,  positive  toward  north 
r>0  ■>  radial  motion  toward  epicenter,  and 

t>0  *>  transverse  motion  to  right  of  propagation  direction. 

Sections  of  Rayleigh  wave  motion  were  then  determined  from  visual  inspection 
of  plots  of  the  radial  and  transverse  components.  Determination  of  the 
particle  motion  of  small  amplitude  high  frequency  motion  in  the  presence 
of  laige  amplitude  low  frequency  motion  is  difficult.  Ideally,  the  radial 
and  transverse  components  should  be  band-pass  filtered  to  separate  the 
frequencies  and  particle  motion  then  determined  for  specified  frequency 
intervals.  However,  due  to  the  computer  time  involved,  this  was  not  feasible 
for  this  study.  Since  only  the  first  higher  Rayleigh  mode  was  present  on 
the  recordings,  the  particle  motion  for  the  frequencies  in1  dved  could  be 
fairly  well  obtained  from  the  unfiltered  radial  and  transverse  components. 


Vertical  component  records  were  convolved  with  101  point,  digital, 
band-pass  filters  described  by  Minear  [1966].  Pass  bands  of  60-100,  25-62, 
16-30,  10-17,  and  7-13  sec  were  used  successively.  The  filtered  data  was 
plotted  by  using  a  Calcomp  plotter.  Period  was  obtained  by  reading  the 
peak-to-peak  period  from  the  Calcomp  plots.  Arrival  time  for  the  period 
was  taken  as  the  time  defined  by 


with  the  quantities  defined  in  the  following  figure. 


tpl  lA  tp2 


Group  velocity  wis  then  obtained  by  dividing  the  epicentral  distance  by 
the  arrival  time. 

Group  velocity  vs.  period  data  for  the  Southern  Appalachians  is 
presented  in  Table  III  and  _  rigs. 8-11. 
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TABLE  111.  RAYLEIGH  WAVE  CROUP  VELOCITIES  FOR  T*IE  SOUTHERN  APPALACHIANS 


EPlCFNTER-J.  ALASKA 

H-m  VTMt->.  ALASKA  STATION-OXFORD,  KISS. 

STATION-Al'I  ANT  A,  GA. 


Fundamental  K.vlelgh  Flrat  Higher  Rayleigh  Fundamental  Rayleigh  Fundament* I  Rayleigh 

iVrlod(aec)  Group  Velocity ( <a/*ec )  Perlodiaec/  Group  Velocity (ka/eec)  Period(eec)  Group  Velocity (ka/eec)  Parlod(aec)  Group  Veloc lty (kn/aec > 


S*» .  on 

3.  72 

8.64 

3.52 

67.04 

3.84 

12.40 

3.01 

52.  75 

1.72 

7.85 

3.56 

54.79 

3.67 

12.40 

2.96 

••1.11 

3.53 

7.85 

3.55 

40.04 

3.37 

11.78 

3.07 

Jb.40 

3.  b2 

7,85 

3.51 

34.54 

3.49 

)  1. 62 

3.05 

3J.H1 

3 . 4  4 

7.07 

3. *9 

31.40 

3.41 

10.52 

2.98 

30.62 

3.38 

7.07 

3.47 

29.05 

3.35 

9.89 

2.97 

26.  b<> 

3.29 

7.07 

3.49 

26.85 

3.14 

9.89 

2.92 

27.  uO 

3.  24 

6.28 

3.46 

24.81 

3.29 

9.89 

2.W 

27.  IK) 

3.19 

6.28 

3.45 

24.65 

3.24 

9.73 

2.95 

24.02 

3.14 

6.91 

3.50 

23.71 

3.09 

9.73 

2.88 

H.31 

3.05 

6.12 

3.54 

23.08 

3.05 

8.32 

2.93 

18.  68 

3.02 

6.12 

3.47 

21.51 

2.97 

8.16 

2.86 

18.  84 

2.96 

5.65 

3.52 

20.25 

3.01 

8.01 

2.89 

18.21 

2.99 

5.34 

3.43 

19.94 

2.94 

18.0b 

2.86 

4.55 

2.53 

19.94 

2.93 

18  Ob 

2.31 

11.37 

2.91 

17.90 

2.91 

18.06 

2.85 

lb. 80 

2.94 

17.74 

2.90 

lb. 17 

2.88 

17.43 

2.88 

IS.  70 

2.83 

16.64 

2.73 

1  r> .  1  7 

3.04 

16  17 

2.83 

14.13 

2.98 

16.17 

2.76 

13.25 

3.02 

16.01 

2.78 

13.03 

3.07 

15.86 

2.80 

12.09 

3.13 

14.97 

2.88 

12.09 

3.00 

14.86 

2.86 

10.99 

3.11 

14.70 

2.83 

1*  99 

3.09 

13.56 

2.81 

9.26 

3.01 

13.19 

3.03 

7.54 

3.01 

12.87 

2.98 

7  44 

i.on 

?  .22 

2.99 

6.59 

2.98 

EPICENTER- JALISCO,  MEX 

EPICENTER -JALISCO,  MEX. 

ST  AT ION- ATLANTA ,  GA. 

STATION- SLACKS 5U1C,  VA. 

Fundamental 

Rayleigh 

Flrat 

Hi,  lyleigh 

Fundai 

■eotel  Rayleigh 

rir.t  Higher  HjyWLh 

Periodic)  Group  Velocity (ka/eec) 

Per lod (etc) 

Group  Velocity (ka/eec) 

Period (aac) 

Croup  Velocity (ka/aac) 

Period(aec) 

Group  Velocity (kn/eec) 

55.  H9 

3.76 

10.05 

3.25 

65.16 

4.11 

18.37 

?.  76 

44.90 

3.52 

9.89 

3.29 

62.02 

3.81 

16.17 

3.84 

18,7? 

3.34 

9.89 

3.22 

50.08 

3.55 

16.01 

3.98 

33.76 

3.20 

8.79 

3.16 

39.41 

3.38 

14.13 

3.32 

27.48 

3.20 

8.32 

2.98 

32.34 

3.7  6 

13.97 

3.23 

22.45 

3.11 

8.01 

3.13 

31.87 

3.26 

13.82 

3.91 

19.94 

3.03 

8.01 

3.05 

29.67 

3.07 

13.66 

3.59 

17.11 

2.97 

8.01 

3.00 

26.38 

2.99 

i3.35 

3.70 

lb. 80 

2.92 

7.85 

3.08 

25.90 

2.92 

13.35 

3.54 

12. 40 

2.89 

7.07 

3.18 

25.40 

2.87 

13.03 

3.49 

12.09 

2.70 

6.28 

3.10 

19.80 

2.77 

12.25 

3.64 

11.93 

2.66 

5.81 

3.03 

19.50 

2.7V 

12.09 

3. 3h 

11.78 

2.85 

18.10 

2.69 

11.93 

3.45 

11.78 

2.82 

14.29 

2.77 

11.93 

3.27 

11.62 

2.78 

13.50 

2.74 

11.78 

3.19 

lo.  99 

2.69 

12.56 

2.70 

11.46 

3  40 

-.0.99 

2.66 

12.09 

2.63 

10.68 

3.13 

10.36 

2.75 

11.93 

2.65 

10.05 

3.1b 

10.  21 

2.75 

11.78 

2.60 

9.73 

3.04 

10.05 

7.72 

9.89 

2.58 

9.26 

3.  In 

10.00 

2.69 

8.01 

2.56 

8.16 

3.03 

9.89 

2.73 

6.  12 

2.55 

8.01 

3.07 

6.90 

2.62 

7.85 

3.01 

H.R«» 

2.6' 

5.97 

.03 

8.20 

2.64 

8.1'#  2.*0 

6.16  2.58 

6.1  6  2.*)* 

7.69  2.*>k 


T  A  ill  .F  IH  (I'OHT’n) 


l 


kph:kktfr-s.  ai,aska 

STATION- BLACKSBURG,  VA. 


Fundamental  Rayleigh 

ttret 

Higher  Rayleigh 

P»rtod(*ec) 

Group  Velocity (k«/»#c) 

Perlod(Mc) 

Group  Vel jcity (km/ sec) 

hh.88 

4.06 

33.13 

4.27 

62.64 

3.88 

32.97 

4.16 

58.88 

3.72 

31.09 

4.43 

49.14 

3.73 

30.93 

4.32 

39.41 

3.62 

30.46 

4.54 

39.41 

3.53 

28.73 

4.09 

31.40 

J  45 

20.72 

4.23 

27.00 

3.20 

20.. 10 

4.35 

25.43 

3.15 

17.11 

4.18 

25.28 

3.10 

15.54 

4.29 

21.82 

3.06 

15.23 

3.94 

19.94 

3.03 

14.44 

3.98 

19.94 

2.97 

13.50 

4.12 

18.68 

2.89 

13.03 

4.20 

18.06 

3.00 

13.03 

4.16 

18.06 

2.94 

13.03 

4.02 

17.58 

2.87 

12.72 

3.90 

17.58 

2.84 

12.72 

3.80 

16.17 

2.89 

12.40 

3.86 

15.86 

2.92 

12.25 

3.83 

15.39 

2.87 

11.93 

3.76 

13.97 

2.84 

11.46 

3.73 

11.62 

3.08 

11.78 

3.70 

11.62 

3.04 

9.69 

3.49 

11.30 

3.06 

8.01 

3.52 

10.21 

3.10 

7.85 

3.46 

9.42 

3.02 

7.69 

3.60 

9.26 

3.00 

6.12 

3.58 

8.32 

2.98 

6.12 

3.54 

8.16 

2,92 

6.12 

3.51 

8.01 

2.93 

6.12 

3.47 

7.85 

3.06 

5.97 

3.55 

7.85 

2.94 

5.65 

3.57 

7.22 

3.00 

7.22 

2.99 

7.22 

2.97 

6.91 

3.07 

6.91 

3.04 

6.75 

3.05 

6.75 

3.03 

6.28 

2.96 

6.28 

2.90 

6.12 

3.02 

5.65 

2.91 

fc- 


5.0  THEORETICAL  RAYLEIGH  DISPERSION  CURVES 

The  basic  Thomson-Haskell  matrix  method  was  used  to  compute 
phase  and  group  velocity  vs.  period  curves  for  layered  earth  models.  A 
computer  program  was  written  to  compute  Rayleigh  wave  dispersion  curves 
and  mode  shape  using  the  modified  formulation  of  the  Thomson-Haskell 
method  presented  by  Dunkin  [1965].  A  program  was  also  written  to  compute 
mode  shape  using  the  Thomson-Haskell  method  A  discussion  and  comparison 
of  the  computation  methods  used  is  given  in  this  section.  Appendix  I  and  II 
contain  detailed  descriptions  of  the  actual  mechanics  of  computation 
and  computer  programming.  Fortran  listings  of  the  programs  are  given  in 
Appendix  III. 

5.1  Thomson-Haskell  Matrix  Method 

As  is  well  known,  the  Thomson-Haskell  matrix  method  consists 
of  evaluating  the  roots  of  a  determinant  formed  by  the  repeated  multiplication 
of  4  x  4  layer  matrices  which  are  functions  of  the  layer  parameters  of  density, 
thickness,  compressional  velocity,  shear  velocity,  as  well  as  phase  velocity 
and  period. 


Using  Haskell's  notation,  the  displacement-stress  matrices  at  the  top 
and  bottom  of  the  m^1  layer  are  given  by 

\ 


where  -  D  F  is  the  m-*-  layer  matrix.  By  repeated  application  of 
(5.1-1),  Haskell  shows  that,  assuming  no  stresses  at  the  free  surface 
.  «  T  =  0)  and  no  sources  at  infinity, 

V*  O 


r  1 

»  • 

F 

1 

!  A  * 

• 

U 

0 

n 

c 

A  ■' 
n 

-  J 

• 

w 

0 

c 

to  * 
n 

0 

U  ' 

n 

0 

L  J 

(5.1-2) 


where  J  «  F  1  a  t  . . .  a, 
n  n-1  1 

a  ,  eliminating  A  '  and  u  ’ 
m  °  n  n 


is  the  matrix  product  of  the  4x4  layer  matrices 
between  the  four  equations  yields 


J22  "  J12  J42  "  J32 

—  ■—  ~  "  -  1  ■  1  s  ■  ■  ■  -  ■  ■■■■ 

Jn  “  J2i  J:i  "  J4i 


(5.1-3) 


Since  the  ,  are  functions  of  phase  velocity  and  wave  number  ,  (5.1-3)  is 

an  implicit  relation  between  c  and  k  and  thus  the  phase  velocity  dispersion 

function.  The  layer  matrix  elements  of  a  are  either  trigonometric  or 

m 

hyperbolic  functions  depending  on  whether  the  phase  velocity  is  greater  than 
or  less  than  the  layer  comp l ess ional  and/or  shear  velocities.  The  multiplica¬ 
tion  of  real  and  imaginary  components  of  matrices  on  a  computer  which  does 
not  have  complex  number  subroutines  would  in  genera.1  add  considerable  complexity 
to  the  problem.  However,  a^  shown  in  Appendix  I,  the  form  of  the  layer 
matrices  leads  to  a  simple  solution  by  which  the  multiplication  of  the  matrices 
with  real  and  imaginary  elements  can  be  accomplished  by  the  multiplication 
of  certain  elements  by  +  1. 
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5.2  Numerical  Difficulties  in  the  Thomson-Haskell  Method 


In  practice,  numerical  computational  difficulties  are  en¬ 
countered  in  the  repeated  matrix  multiplication  required  to  evaluate  the 
roots  of  (5.1-3).  These  difficulties  are  encountered  as  the  product  of 
kH,  where  k  is  the  wave  number,  and  H  is  the  total  thickness  of  the  layered 
earth  model,  becomes  large.  Dorman,  Ewing,  and  Oliver  [I960]  have  used 
an  upper  limit  of  about  30  for  kH.  When  the  value  of  kH  reaches  about  30, 
the  number  of  layers  can  be  reduced  and  the  computation  continued  with  a 
reduced  thickness.  Little  error  is  introduced  by  thir  technique.  However, 
for  higher  modes  and  hence,  higher  frequencies,  the  product  kH  may  be  relatively 
large  even  for  layered  earth  models  of  small  total  t'nicxness,  H. 

Dunkin  [1965]  has  shown  the  numerical  difficulties  are  caused  by  the 
computation  of  large  exponentials  and  a  resulting  loss  of  singificant 
figures.  Dunkin' s  development  is  briefly  repeated  below  in  order  to  show 
the  effect  of  loss  of  significance  due  to  the  addition  of  large  and  small 
quantities  on  the  Haskell  matrix.  The  agrument  is  applied  directly  to  the 
Haskell  dispersion  equation  (5.1-3)  rather  than  to  the  secular  equation  used 
by  Dunkin. 

Let  the  scalar  and  vector  potential  functions  for  an  elastic  body  have 
the  form 


4>n  =  expik(ct-x)  [A^  exp(ikz  exp(-ixz  Vc2/a2-l)  ] 


expik(ct-x)  [4>n+  +  <j»n  ] 


n 


expik(ct-x)  [Cn  exp(ikz  V/c"/B2 -1)  +  exp(-ikz  Vc2/B2-1)  ] 


(5.2-1) 


expik(ct-x)  -  i|>  ] 
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Using  (5.2-1)  and  the  equation 


S  (u,v,w)  «  ^  *  <f>  +  $  x  (ip1,  ^3) 


the  displacement-stress  vector, 


S^,  can  be  expressed  as 


S  (2)  -  T  *  (2) 

n  n  n 


where 


Vz) 


r 

i 


and  is  a  4  x  4  matrix  function  of  c,  k,  and  the  layer  parameters. 


Taking  the  origin  at  the  2^  ^  interface  (5.2-1)  gives 


n-1 


A  +  B 


n  n 


(zn-l) 


+  ^n 


) 


n-1 


C  +  D 
n  n 


>  +  *n 


At  the  2  interface 
n 


(5.2-2) 


(5.2-3) 
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4>  -  A  expikd  r  +  B  exp-ikd  r  =  i  (z  )  +  <k  (z  ) 

n  n  n  an  n  n  an  an  nn 

ip  =  C  expikd  r„  +  D  exp-ikd  r„  =  ip  +  (z  )  +  'i'  (z  ) 

n  n  r  n  8n  n  ngn  nn  nn 


The  relation  betv/een  <t>  (z  )  at  the  z  interface  and  4>  (z  1)  at  the  z  .. 

n  n  n  n  n-1  n-1 


interface  is  then 


4> 

n 


<-„> 


E  *  (z  ,) 
n  n  n-1 


(5.2-4) 


where 


E 

n 


expikd  r 
r  n  an 

0 

0 

0 

0 

explkV8n 

0 

0 

(5.2-5) 

0 

0 

exp-ikd  r 
r  n  an 

0 

0 

0 

0 

exp-ikd  rQ 
r  n  Bn 

At  the  n-1  interface 


S 

n 


(zn-l> 


T  4>  (z  ,) 
n  n  n-1 


or 


4> 

n 


) 


Sn(Zn  1> 

n  n-i 


(5.2-6) 


Now,  by  the  boundary  conditions  of  the  continuity  of  stress  and  displacement 


S 

n 


(zn> 


T  t  (z  ) 
n  n  n 


(5.2-7) 


Substituting  (5.2-4)  for  4>  (z  )  and  (5.2-6)  for  4>  (z  -.)  in  (5.2-7)  gives 

n  n  n  n-i 
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-1 


(5.2-8) 


S  (z  )  =  T  E  T  S  (2 

n+1  n  n  n  n  n  n-1 


This  equation  is  equivalent  to  Haskell's  equation 


S  (z  )  ■  a  S  (z  -)  , 

n+1  n  n  n  n-1  ’ 


(5.2-9) 


with  a 


n 


T  E  T 
n  n  n 


-1 


By  (5.2-8)  the  displacement-stress  vector  is  converted  into  $  ,  continued 

n 

through  the  layer  z  by  E  ,  and  ccnverted  back  into  S  (z  )  at  the  interface 
0  1  n  J  n  n  n 

n+1  which  is  equal  to  Sn+^  (z^).  The  Haskell  layer  matrix  carries  the  dis¬ 
placement-stress  vector  from  the  n^  interface,  through  the  layer,  and  across 
the  n+1  interface  in  one  operation.  Eq.  (5.2-8)  brings  into  evidence  the 
effect  of  the  "continuing"  matrix  En« 

Consider  the  matrix  linking  the  displacement-stress  vectors  at  the  free 
surface  and  the  last  layer  of  an  assumed  layered  sequence. 


S 


n-1 


n-1 


G.S  =  PS 
1  o  o 


(5.2-10) 


where 


G 

n 


T  E  T 
n  n  n 


In  the  Haskell  formulation,  the  matrix  from  which  the  dispersion  relation 

is  obtained  is  given  by  J  =  F  and  in  the  Dunkin  formulation  this  matrix 

n 

is  T  ^P.  However,  considerations  of  the  numerical  evaluation  of  the  P 
n 

matrix  will  yield  results  valid  to  both  developments  since  the  T^  ^  or 
^  do  not  contain  exponential  powers. 
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Let  P  be  written  as 


P  -  A  T  F.  T  B 
in  m  m 


(5.2-11) 


where 


A  =  G  ,  . . .  G  . 

n-x  m+1 


B  =  G  ....  G. 

m-1  1 


Using  the  definitions  of  T  and  E  ,  it  can  be  shown  that  the  components 

mm 

of  P  are  of  the  form 


P..  *  B  expikd  r  +  C  expikd  rQ  +  D.  .  exp^ikd  r  +  E..  exp^ikd  r 

^  -j  ^  -j  r  m  rvm  1  1  1  m  ftm  i  r  m  n»m  11  r  m 


m  8m 


m  am  ij 


m  8m 


(5.2-12) 


For  the  Haskell  development,  is  of  the  form 


F11  0  F23  0 

0  F22  0  F24 

F31  0  F33  0 

0  F42  0  F44 


which  gives  for  the  two  components  and  ^22 


J12  "  F11P12  +  F13P32  *  and 


(5.2-13) 


J22  "  F22P22  +  F24P42 
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tn 

Now  suppose  that  for  the  m*“  layer  r  and  rQ  are  negative  imaginary 

am  pm 

(c<a)  so  tha  expidk  r  and  expikd  rfl  may  be  large  depending  on  the  value 

m  am  m  pm 

of  k.  If  the  exponential  term  is  large  enough,  the  effect  of  the  smaller 
terms  in  will  be  neglected  in  computing  the  P  ,  because  of  loss  of 
significance.  In  the  evaluation  of  the  roots  of  (5.1-3)  the  difference 
(Jl2  -  must  be  taken.  Although  the  P  are  large,  their  differences 

may  be  small.  Therefore,  terms  which  were  lost  because  of  loss  of  significance 
in  computing  the  P^  would  be  important  in  the  difference  of  and 
Mode  shape  is  computed  from  repeated  applications  of  (5.1-1)  using  the 
starting  values  of  u^  and  *wq  from  (5.1-3).  Therefore,  the  same  problem 
with  loss  of  significance  is  inherent  in  the  Haskell  method  of  computing  modal 
shape. 

5.3  Dunkin  Modification  of  the  Haskell  Method 

Dunkin  derives  the  secular  or  period  equation  in  the  form 


Det 

*u  ■  0 

(5.3  1) 

Where 

R11 

R  *  R 

R21 

R12 

R22 

=  T  "1  G  .  ... 
P  P-1 

G1 

(5.3-2) 

He  shows  that  Det  can  be  expanded  as  a  product 

of  the  second  order 

subdeterminants  of  T  ^ 

P 

and  G 

P 

yielding 

Det  = 

tp-i 

(12  p-1  jaK  1 

|ab  8  |cd  *'*  8 

|ef 

1 1 2  ’ 

(5.3-3) 

25 


1) 


where  g  ^  is  the  second  order  subdeterminant  of  G  involving  rows 

Ikf  P 

i  and  j  and  columns  k  and  JL.  Dunkin  has  shown  that  by  using  algebraic 
expressions  for  the  subdeterminants  of  the  ,  numerical  difficulties 
with  loss  of  significance  can  be  avoided  since  the  products  of  like  ex¬ 
ponentials  normally  occurring  in  the  secular  function  are  excluded  at  the 
start.  Products  of  unlike  exponentials  for  a  given  layer  effectively  in¬ 
crease  the  magnitude  of  Det  R^.  To  prevent  machine  overflow,  the  secular 
function  can  be  divided  by  the  two  largest  exponents  when  these  exponenets 
become  real  and  the  exponential  expression  becomes  hyperbolic.  This  results 
in  no  loss  of  significance. 

Explicit  expressions  for  the  gjj^  and  the  g  are  given  in  Appendix 
II  for  real  frequencies  and  wave  numbers.  These  are  slightly  different 
from  the  definitions  of  Dunkin, since  he  assumes  complex  frequencies. 

Mode  shapes  are  computed  using  the  following  relation  of  Dunkin  dis¬ 
cussed  in  Appendix  II  , 
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5.4  Computational  Procedure 

The  Dunkin  method  was  programmed  in  Fortran  II  for  the  Bunker- 
Ramo  340  computer.  Equation  (5.3-1)  was  used  for  the  determination  of  Ray¬ 
leigh  wave  dispersion  curves.  Equation  (5.3-4)  was  used  to  compute  the 
mode  shape  once  the  roots  of  (5.3-1)  were  obtained.  A  double  precision 
program  was  written  in  Fortran  II  for  the  IBM  360-75  using  the  Haskell  method 
for  computing  mode  shape.  Equation  (5.1-1)  was  used  for  this  computation. 
Figs.  5  and  6  are  flow  charts  of  the  computer  program  FLATRAY  used  in  the 
computation  of  Rayleigh  dispersion  and  modal  shape.  A  computer  listing  of 
the  program  is  given  in  Appendix  III. 
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Fig.  5.  Flow  chart  of  FLATRAY 
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The  technique  used  for  finding  the  roots  of  (5.3-1)  consists  of 
two  steps.  In  order  to  determine  roots  for  different  modes,  a  gi.’en  phase 
velocity,  cq,  and  starting  wave  number,  k^Q  ,  are  specified.  is  then 

incremented  using  the  constant  cq  value  until  a  root  is  bracketed.  A  two 
point  interpolation  scheme  is  then  used  until  the  difference  between  two 
values  of  k  which  successively  bracket  the  root  is  less  than  an  input  value. 

If  more  than  one  mode  is  to  be  investigated,  k  is  incremented  from  its 
value  at  the  last  root  found  using  the  same  value  of  cq  until  the  next  root 
is  found.  Thus,  the  roots  along  a  constant  c  curve  are  found  which  correspond 
to  different  modes.  In  order  to  define  a  particular  mode,  computation  begi  ,s 
at  cq  and  the  k^Q  corresponding  to  the  desired  mode.  C  is  decremented  and 
k  incremented  by  values  specified  as  input  parameters.  K  is  varied  at  the 
new  -  Ac  until  the  root  is  bracketed.  Two  point  interpolation  is  then 
used  until  k  is  obtained  with  the  desired  accuracy.  The  process  of  de¬ 
crementing  c  and  incrementing  k  is  continued  until  three  points  on  a  given 
modal  dispersion  curve  are  found.  A  three  point  Gregory-Newton  interpolation 
scheme  is  then  used  to  estimate  the  next  root  on  the  curve.  The  process  of 
interpolation  and  bracketing  continues  until  a  dispersion  curve  is  defined 
to  some  minimum  specified  value  of  c. 

Group  velocity  values  are  computed  by  perturbing  c  a  small  amount 
from  a  value  at  which  a  corresponding  k.  udS  been  found.  K  is  found  for 
the  perturbed  c  and  a  two  point  difference  scheme  used  to  evaluate  group 
velocity,  U  ,  according  to 
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From  Table  VI  i  t  can  be  seen  that  the  group  velocity  curves  computed  by 
the  Dunkin  method  for  many  layered  models  agree  with  those  computed  for 
models  with  reduced  thickness  to  within  a  few  te:  ths  of  a  per  cent. 

However,  the  section  must  be  taken  thick  enough  to  include  the  entire 
depth  to  which  appreciable  particle  motion  extends.  This  is  illustrated 
by  the  first  higher  mode  for  the  Gutenberg-Birch  II  model.  The  period 
value  at  5.00  km/sec  is  50.0301  sec  for  the  400  km  section,  and  46.4060 
sec  for  the  1000  km  section.  However,  reference  to  the  displacement  vs. 
depth  curves  in  Appendix  IV  shows  that  the  ^0  km  section  does  not  include 
the  total  depth  to  which  vertical  particle  motion  extends  at  this  period. 

Although  thick  sequences  can  be  used  to  compute  shorter  period  group 
velocity  curves,  it  is  considerably  more  economical  in  computer  time  to 
use  thinner  sequences  of  fewer  layers.  Displacement  depth  curves,  such 
as  those  in  Appendix  IV,  can  be  used  to  indicate  the  necessary  total 
thickness  to  be  used  at  given  periods. 

Modal  Shape 

After  a  point  (c,k)  was  found  on  a  given  dispersion  curve,  the 
mode  shape  for  the  given  (c,k)  was  computed  using  (5.3-4).  Horizontal 
and  vertical  displacements  vs.  depth  values  were  then  punched  out  on  cards 
to  be  used  in  computing  the  variation  of  phase  velocity  due  to  variation 
in  layer  parameters  from  INTEGRAL.  Modal  shape  was  computed  using  both 
the  Dunkin  and  Haskell  methods.  A  double  precision  program  was  used  in 
computing  modal  shape  by  the  Haskell  Method.  Table  IV  shows  the  comparative 
results  for  the  two  methods.  Values  for  horizontal  and  vertical  particle 
amplitudes  agree  very  well  ( < . 03%)  for  the  first  four  layers  in  each  case 
considered  in  Table  IV.  After  this,  the  differences  between  the  two  methods 
increase  rapidly.  The  rapid  increase  of  particle  amplitude  with  depth  in 
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TABLE  IV.  Fundamental  Rayleigh  mode  vertical  and  horizontal 
Particle  amplitudes  computed  by  Dunkin  and  Haskell 
Methods  for  the  Gutenberg-Birch  II  model.  Displacement 
normalized  to  the  vertical  displacement  at  the  surface 


Period  =  25.  ,440 
Phase  Velocity  *  3.8000 
7  Layer  model 


HASKEL 

Vertical 

Section  Thickness  * 

L  METHOD 

Horizontal 

140  km 

DUNK 

Vertical 

1.000000 

.6950871 

1.000000 

.886979 

.0084581 

.886972 

.518503 

.0569591 

. 518524 

.303885 

.1063181 

.304132 

.164817 

.0761211 

.062881 

.078571 

.0509691 

.014165 

.018264 

.3484311 

.002693 

Period  -  25.1380 

Phase  Velocitv  3.8000 
17  Layer  mdel 

1.000000 

Section  Thickness 

-.6950441 

=  400  km 

1.000000 

.886837 

.0086111 

.886883 

.518193 

.0569161 

.518278 

.303773 

.1056621 

.303749 

.165958 

.0734501 

.062646 

.084840 

.0412981 

.013984 

.043072 

.0154241 

.002541 

.034665 

-.0212121 

-.000083 

.087474 

-.1191061 

-.000634 

.324986 

-.4357481 

-.000713 

1.117859 

-1.4545031 

-.000685 

4.078959 

-4.9270861 

-.000643 

13.757256 

-16.1905901 

-.00059* 

45. >49210 

-53.0181841 

-.000553 

150.983339 

-173.2435981 

-.C30512 

Horizontal 


-.6950871 

.0084451 

.0568251 

.0400851 

.0126391 

.0034631 

.0009951 


-.6950441 
.0085831 
.0569481 
.(401191 
.0126671 
.0034971 
.0010301 
.0004061 
.0002461 
.0001991 
.0001791 
.0001651 
.0001531 
. 00U1431 
.0001331 
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table  IV  (Cont'd) 


Period  =  25.1380  (Cont'd) 
Phase  Velcoity=3. 8000 
17  Layer  model 


Section  Thickness  =  403  km 


HASKELL  METHOD 

1 1  DUNKIN.  METHOD 

Vertical 

Horizontal 

Vertical 

Horizontal 

8.176407 

-567.5807251 

(  -.000474 

.0001231 

10070.040164 

l  -11362.1297281 

1  -.001556 

.0000261 

Period  =  106.4400 

Phase  Velocity=4. 2000 

35  Layer  model 

Section  Thickness 

=  2898  km 

1.000000 

-.8380621  | 

1.000000 

-.8380621 

1.047271 

j 

-.5990881  ‘ 

i 

1.047270 

-.5990121 

1.054375 

-.4060961  | 

j 

1.054340 

-.4059721 

1.048108 

-.1909051 

1.048030 

-.1826711 

1.01G899 

-.0527151 

.968018 

-.0486231 

.954535 

.0438221 

.883244 

.0393891 

.886795 

.1098891 

.793526 

.0953051 

. 81252C 

.1520151 

.703175 

.1274421 

.735438 

.1759571 

.615440 

.1425871 

.658577 

.1869151 

.532770 

.1463311 

.583974 

.1879691 

.456412 

j 

.1420871 

.513373 

.1825901 

.387350 

.1331501 

.447853 

.1730011  1 

j 

.325880 

.1216331 

.387824 

.1607891  j* 

.271911 

.1088621 

.333572 

.1476171  !i 

|  • 

.225107 

.0961331 

.284944 

- 1341261 

.184894 

.0839571 

. 184929 

.0991591 

. 115100 

.0486511 

• 1’ 4067 

.0723331 

.055227 

.0274971 

.024591 

.0434351 

.008158 

.0065751 

-.005338 

.0396271 

-.001496 

.0047631 

-.034168 

.0439991 

j 

-.004566 

.0042311 

-.117644 

.0864391  | 

j 

-.011717 

.0032201 

-.299466 

.1955751  ] 

-.010744 

.002691-l 
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TABLE  IV  (Cont'd) 


Period  ■  106.4400  (Cont'd) 
Phase  Velocity»4. 2000 
35  Layer  model 


Section  Thickness  ■  2898  km 


HASKELL  METHOD 


Vertical 

Horizontal 

-.734195 

.4387981 

-1. 73:194 

.8698931 

-7.174500 

-1.4270811 

67.316536 

-162.8932311 

2727.220784 

-3866.6973961 

58523.227239 

-V2877. 2524311 

.108  x  10? 

-.127  x  107  1 

.188  x  108 

-.214  x  108  1 

.319  x  109 

-.355  x  109  1 

.525  x  1010 

-.576  x  1010i 

.855  x  1011 

-.928  x  10Ui 

.329  x  1012 

-.352  x  1012i 

DUNKIN  METHOD 


Vertical 

Horizontal 

-.009416 

.0022971 

-.008134 

.0019341 

-.023747 

.0005461 

-.006543 

.0001411 

-.001763 

. 0000361 

-.000463 

.0000091 

-.000119 

.0000021 

-.000030 

.0000001 

-.000007 

. OOOOOOi 

-.000002 

. OOOOOOi 

-.000000 

. OOOOOOi 

-.000000 

.0000001 
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the  Haskell  Method  results  from  the  repeated  multiplication  of  4  x  4 
layer  matrices  and  gradual  loss  of  significance  in  the  matrix  multi,  ' ication. 
Since  the  Haskell  Method  starts  at  the  top  layer  and  works  down, and  the 
Cunkin  Method  starts  at  the  bottom  layer  and  works  up,  the  agreement  of 
the  two  methods  in  the  near  surface  layers  indicates  that  the  Dunkin  Method 
is  yielding  correct  displacement  values  over  the  entire  layered  sequence. 

In  some  cases  when  relalive  high  frequency  points  on  a  dispersion  curve 
were  being  computed  with  a  many  layered  model,  displacements  computed  by 
the  Dunkin  Method  showed  some  slight  perturbations  with  depth  rather  than 
a.  smooth  decrease.  Displacement  values  also  tended  to  change  signs  as 
they  decreased  to  very  small  quantities  with  depth.  This  is  seen  in  the 
vertical  displacements  for  the  35  layer  case  in  Table  IV.  Horizontal  displace¬ 
ment  curves  generally  have  one  more  lobe  than  the  corresponding  vertical  dis¬ 


placement  curves. 

Modal  shapes  for  the  Gutenberg-Birch  model  are  shown  graphically  in 
Appendix  IV. 

Earth  Flat tening  Approximation 


The  earth  flattening  approximation  introduced  by  Alterman,  Jarosch, 
and  Pekeris  [1961]  was  used  to  modify  the  layer  velocities;.  As  has  been 
shown  by  Kovach  &  Anderson  [1964], the  effect  of  sphericity  is  not  negligible 
even  for  higher  modes.  The  linear  increase  in  velocity  introduced  in  the 
earth  flattening  approximation  is  specified  by  the  parameter 


The  value  of  r  for  a  layer  was  taken  as  the  radius  to  the  center  of  the 
m  J 

layer;  a  is  the  mean  radius  of  the  earth,  6371  km.  Layer  velocities  approxi¬ 
mately  corrected  for  sphericity  are  then  given  by 
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5.5  Variation  of  Phase  Velocity  with  Layer  Parameters 

The  energy  integrals  for  elastic  wave  propagation  [Meissner, 
1926;  Jeffreys,  1934]  have  been  used  by  several  authors,  notably  Anderson 
[1964]  and  Takeuchi  and  Dorman  [1964]  to  derive  explicit  relations  between 
the  variation  of  phase  velocity  and  the  variation  of  layer  parameters. 

Necessary  data  for  the  evaluation  of  the  partial  derivatives  or  phase  velocity 
with  respect  to  layer  parameters  are  horizontal  and  vertical  particle  amplitude 
vs.  depth  values. 

For  Rayleigh  waves,  the  potential  and  kinetic  energy  averaged  over  a 
cycle  are 


4T 


2  2  2 
pm  (u  4  w  )  dz 


4V 


[A(!  u-w')2  4  y(2k2u2  +  2w'2)  4  k2w2  4  u'  2  +  2ku'w]  dz 


where 

u  *  horizontal  displacement 

w  =  vertical  displacement 

u1  ■  ~  ,  and  the  integration  extends  over  the  entire  depth. 

(7  Z 

Using  the  fact  that  the  kinetic  and  potential  energy  averaged  over  a  period 
are  equal  we  obtain 


35 


(5.5-1) 


where 
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For  a  layered  sequence  of  n  layers,  one  can  define 
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and  similarly  for  the  other  integrals.  Thus, 
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A  perturbation  of  a  layer  parameter  in  the  m —  layer  will  cause  a 
in  the  integral,  1^  ,  for  the  entire  layered  sequence  of 
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Differentiating  (5.5-1)  with  respect  to  the  layer  parameters,  the  partial 
derivatives  of  c  with  respect  to  the  layer  parameters  are  obtained.  Thus, 
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where 

D  -  k[2k(I2  +  I5)  +  I3  +  I6] 

with  the  integration  of  the  integrals  in  D  extending  over  the  entire  layered 
sequence. 
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Assuming  p,a,B  as  independent  and  p,A,y  as  dependent 
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gives 
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The  group  velocity  can  be  expressed  in  terms  of  the  integrals  1^  by 
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Equations  (5.5-6),  (5.5-8),  (5.5-9),  (5.5-10),  (5.5-11),  (5.5-12), 
and  (5.5-13)  with  the  definitions  of  (5.5-2)  were  programmed  for  the  IBM 
360-75  computer.  The  integrals  of (5.5-2)  were  evaluated  numerically  using 
polynomial  approximations  to  the  particle  displacements  obtained  from  the 
modal  shape  calculations.  A  sliding  fitting  procedure  was  used  in  determining 
the  polynomials.  In  thir  procedure,  a  polynomial  is  fitted  to  say  n  points 
at  the  depths  z^,  ,  ...  ^i+n  *  and  the  integrals  and  their  derivatives 

evaluated  over  the  interval  z^+n  “  The  polynomial  fit  is  then  shifted 

to  drop  one  layer  and  pick  up  one  layer,  i.e.,  to  the  points  at  depths 
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,  z^jr2  *  ***  zi+l+n  ’  By  using  polynomial  approximations,  the  integrals 
of  the  derivatives  are  simply  the  integrals  of  the  derivatives  of  the 
polynomials.  Third  degree  polynomials  were  found  to  give  adequate  fits  to 
the  displacement  data. 

A  Fortran  listing  of  the  program,  INTEGRAL,  for  computing  the  partial 
derivatives  of  phase  velocity  ,c  ,  with  respect  to  layer  parameters  i9 
given  in  Appendix  III,  Results  of  computation  of  the  partial  derivatives 
for  the  Gutenberg-Birch  IX  continental  model  are  given  in  Appendix  IV. 
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6.0  CRUSTAL  AND  UPPER  MANTLE  STRUCTURE  IN  THE  SOUTHEASTERN  UNITED  STATES 


Gravity,  local  travel-time,  and  Rayleigh  wave  dispersion  data 
were  used  to  determine  crustal  and  upper  mantle  structure  for  the  South¬ 
eastern  United  States.  Results  of  the  Rayleigh  dispersion  study  and  a 
comparison  with  the  gravity  and  travel-time  results  are  given  in  this 
section. 

Observed  dispersion  curves  could  be  constructed  only  for  periods  less 
than  about  50  seconds.  Therefore,  crustal  structure  was  concentrated  on 
and  the  upper  mantle  structure  below  70  km  was  assumed  to  be  that  of  the 
Gutenberg-Birch  II  continental  model.  The  Gutenberg-Birch  II  crustal 
structure  was  used  as  the  basic  model  for  estimating  the  Southern  Appalachian 
structure.  Variations  of  phase  velocity  with  layer  parameters  computed  for 
the  Gutenberg-Birch  II  model  (See  Appendix  IV)  were  used  to  vary  this  basic 
model  to  yield  dispersion  curves  fitting  the  observed  data.  Velocity  and 
density  structure  of  the  models  considered  are  given  in  Thble  V  and  Fig. 

7.  Values  of  the  "earth  flattening"  velocities  for  the  Gutenberg-Birch  II 
are  also  given  to  indicate  the  effective  increase  of  velocity  with  depth. 

Fundamental  and  first  higher  Rayleigh  mode  group  velocity  vs.  period 
data  observed  in  the  Southern  Appalachians  (Table  III)  are  shown  in  Figs. 
8-11.  Rayleigh  wave  dispersion  curves  computed  as  described  in  Section 
5.0  are  given  in  Table  VI  and  Figs.  8-11  for  the  models  Gutenberg-Birch  II, 
310,  314,  315,  and  320  defined  in  Table  V. 

Group  velocity  curves  for  waves  traveling  approximately  perpendicular 
(perpendicular  waves)  to  the  Appalachians  are  quite  similar  (Figs.  8-11). 

They  all  have  a  local  minimum  of  about  2.85  km/sec  at  a  period  of  17  sec, 
and  a  local  maximum  of  about  3.10  km/sec  at  around  12  sec.  For  periods 
shorter  than  about  10  seconds  the  curves  flatten  out  at  about  3.0  km/sec. 
First  higher  Rayleigh  mode  curves  indicate  a  broad  minimum  of  3.5  km/sec 
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TABLE  V. 


Layer  Parameters  for  Crustal  and  Upper  Mantle  Models 


GUTENBERG-BIRCH  II 


FLAT  EARTH 

EARTH  FLATTENING 

6  p 

d 

a 

8  p 

a 

6.1400 

6.5800 

8.0800 

7.8700 

7.8000 

7.8300 

7.8900 

7.9400 

8.0000 

8.0600 

8.1200 

8.2000 

8.2700 

8.3500 

8.4300 

8.5100 

8.7500 

9.0000 

9.4900 

9.7400 

9.9900 

10.5000 

10.9000 

11. 3000 

il.4000 

11.8000 
12.0500 

12.3000 

12.5500 

12.8000 
13.0000 

13.2000 
13.4500 

13.7000 

13.6500 


3.5500 

3.8000 

4.6000 

4.5100 
4.4500 

4.4200 

4 . 4000 
4.3900 

4.4000 

4.4200 
4.450" 
4.480C 
4.5200 

4.5700 

4.6100 
4.6600 
4.8100 

4 . 9500 
5.2200 
5.3600 

5.5000 

5.7700 
6.0400 

6.3000 
6.3500 

6.5000 

6.6000 

6.7500 

6.8500 

6.9500 
7.0000 
7.1000 

7 . 2000 

7 . 2500 

7 . 2000 


2.7500 

2.9000 

3.5700 

3.5100 

3.4900 

3.5000 

3.5100 
3.5300 

3.5500 
3.5600 

3.5800 

3.6100 
3.6300 

3.6500 
3.6800 

3.7000 

3.7700 

3.8500 
4.0000 
4.0700 
4.1500 

4.3000 

4.4200 
4.5400 

4.5700 
4.6900 

4.7700 

4.8500 
4.9200 
5.0000 
5.0600 
5.1200 
5.1900 

5.2700 

5.2500 


19.00 
19. 0C 
22.00 
20.00 
2C.00 
20.00 
20.00 
20.00 
20.00 
20.00 
20.00 
20.00 
20.00 
20.00 
20.00 
50.00 
50.00 
100.00 
50.00 
50.00 
100.00 
100.00 
100,00 
100.00 
200.00 
200.00 
200.00 
200.00 
200.00 
200.00 
200.00 
200.00 
200.00 
98.00 
60 


6.1486 

6.6090 

8.1422 

7.9574 

7.9115 

7.9670 
8.0541 
8.1314 
8.2192 
8.3074 
8.3961 
8.5067 
8.6074 
8.7191 
8.8313 

8.9670 
9.2969 
9.6840 

10.3422 

10.7062 

11.1249 

11.9007 

12.5775 

13.2798 

13.7780 

14.8243 

15.7602 

16.7759 

17.8824 

19.0925 

20.3437 

21.7245 

23.3411 

24.7819 

25.0395 


3.5550 

3.8167 

4.6354 

4.5601 

4.5136 

4.4974 

4.4915 

4.4958 

4.5206 

4.5557 

4.6013 

4.6476 

4.7044 

4.7720 

4.8294 

4.9102 

5.1106 

5.3262 

5.6888 

5.8917 

6.1248 

6.5397 

6.9696 

7.4038 

7.6746 

8.1660 

8.6321 

9.2063 

9.7606 

10.3666 

10.9543 

11.6852 

12.4949 

13.1145 

13.2077 


2.7500 

2.9000 

3.5700 

3.5100 

3.4900 

3. 5000 

3.5100 
3.5300 

3.5500 
3.5600 

3.5800 

3.6100 
3.6300 

3.6500 
3.6800 

3.7000 

3.7700 

3.8500 
4.0000 
4.0700 
4.1500 

4.3000 

4.4200 
4.5400 

4.5700 
4.6900 

4.7700 

4.8500 
4.9200 
5.0000 
5.0600 
5.1200 
5.1900 

5.2700 

5.2500 


d 

19.00 
19.00 
22.00 
20.00 
20.00 
20.00 
20.00 
20.00 
20.00 
20.00 
20.00 
20.00 
20,00 
20.00 
20.00 
50.00 
50.00 
100.00 
50.00 
50.00 
100.00 
100.00 
100.00 
100.00 
200.00 
200.00 
200.00 
200.00 
200.00 
200.00 
200.00 
200.00 
200.00 
98.  on 
oo 
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TABLE  V.  (Cont'd) 


MODEL 

,  310 

MODEL  314 

a 

3 

k 

d 

a 

3 

P 

d 

5.8800 

3.3800 

2.6700 

10.00 

5.8800 

3.3800 

2.6700 

10.00 

6.1400 

3.5500 

2.7600 

10.00 

5.6000 

3.2400 

2.7600 

10.00 

6.5800 

3.8000 

2.9000 

10.00 

6.1000 

3.5000 

2.9000 

10.00 

6.5800 

3.8000 

2.9000 

10.00 

6.6000 

3.8000 

2.9000 

10.00 

8.0800 

4.6000 

3.5700 

20.00 

7.0000 

4.1000 

3.1000 

10.00 

8.0800 

4.6000 

3.5700 

20.00 

Same  as 

Gutenberg-Birch  II 

to  400  km 

Same  as 

Gutenberg-Birch  II 

to  810  km 

MODEL  315 

MODEL  320 

a 

3 

P 

d 

a 

3 

P 

d 

5.8800 

3.3800 

2.6700 

10.00 

5.8800 

3.3800 

2.6700 

10.00 

6. 1400 

3.5500 

2.7600 

10.00 

5.8600 

3.3800 

2.7500 

10.00 

6.5800 

3.8000 

2.9000 

10.00 

6.1000 

3.5000 

2.9000 

10.00 

6.5800 

3.8000 

2.9000 

10.00 

6.6000 

3.8000 

2.9000 

10.00 

7.0000 

4.1000 

3. 1000 

10.00 

7.0000 

4.1000 

3.1000 

10.00 

8.0800 

4.6000 

3.5700 

20.00 

8.0800 

4.6000 

3.5700 

20.00 

Same  as  Gutenberg-Birch  II  to  810  km  Same  as  Gutenberg-Birch  II  to  810  km 


(WM)  Hld30 


CRUSTM  AND  UPPER  MANTLE  MODELS 


TABLE  VI.  Rayleigh  wave  dispersion  curves  for  Gutenberg-Birch  II  models, 

and  models  310,  314,  and  315 

GUTENBERG-BIRCH  II 


FUNDAMENTAL  RAYLEIGH  MODE 


Section  Thickness  2898.00 


c  k 


4.5000 

.00850 

4.4000 

.00961 

4.3000 

.01131 

4  1999 

.01405 

4.0999 

.01984 

3.9999 

.03815 

3.8999 

.05419 

3.7999 

.06624 

3.699S 

.07815 

3.5998 

.09202 

3.4998 

.10996 

3.3998 

.14161 

T  U 


164.2671 

148.5391 

3.6990 

129.1709 

3.7695 

106.4402 

3.8361 

77.2556 

3.9112 

41.1794 

3.7346 

29.7316 

3.4039 

24.5623 

3.1759 

21.7314 

3.0493 

18.9673 

2.9510 

16.3264 

3.0336 

13.0509 

3.1084 

Section  Thickness  4l)0.00 


4.0000 

.03801 

3.9000 

.05394 

3.8000 

.06578 

3.7000 

.07808 

3.5999 

.09207 

3.4999 

.11010 

3.3999 

.14220 

41.3293 

29.8693 

3.4253 

25.1380 

3.2096 

21.7493 

3.0305 

18.9572 

2.9639 

16.3054 

3.0331 

12.9958 

3.1018 

Section  Thickness  140.00 


4.0000 

.03730 

3.9000 

.05409 

3.8000 

.06376 

3.7000 

.07762 

3.5999 

.09150 

3.4999 

.11001 

3.3999 

.14249 

42.1107 

29.7861 

3.4066 

25.1423 

3.1768 

21.8791 

3.0923 

19.0747 

2.9947 

16.3184 

3.0426 

12.9695 

3.0588 

FIRST  HIGHER  RAYLEIGH  MODE 


T 


143.8940 

123.5369 

99.2251 

65.9224 

37.1515 

28.3452 

24.0856 

21.0107 

18.3398 

15.5185 

11.9219 


28.4013 

24.1929 

21.0535 

18.3140 

15.4998 

11.9019 


28.3880 

24.2579 

21.0892 

18.3862 

15.49n 

11.8907 


Section  Thickness  2898.00 


c 

k 

1 

U 

T 

6.0000 

.01020 

102.6767 

5.9000 

.01091 

97.6155 

4.503C 

96.2532 

5.8Cx) 

.01174 

92.2541 

4.4644 

90.8938 

5.6999 

.01269 

86.8711 

4.4340 

85.5044 

5.5999 

.01377 

81.4673 

4.4137 

80.0789 

5.4999 

.01503 

75.9839 

4.3971 

74.5686 

5.3999 

.01653 

70.4033 

4.3821 

68.9585 

5.2999 

.0183?. 

64.7021 

4.3720 

C3.2178 
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TABLE  VI.  (Cont'd) 


FIRST  HIGHER  RAYLEIGH  MODE  (Cont'd) 


Section  Thickness  2898.00  (Continued) 


c 

k 

T 

U 

5.1998 

.02053 

5*  578 

4.3631 

5.0998 

.02331 

52.^603 

4.3539 

4.9999 

.02690 

46.7207 

4.3448 

4.8999 

.03172 

40.4258 

4.3363 

4.7999 

.03854 

33.9652 

4.3316 

4.6998 

.04900 

27.2834 

4.3340 

4.5998 

.06780 

20.1477 

4.3586 

Section  Thickness  1000.00 

5.4000 

.01661 

70.0684 

5.3000 

.01837 

64.5405 

4.4587 

5.2000 

.02062 

58.5864 

4.3608 

5.0999 

.02345 

52.5303 

4.3509 

4.9999 

.02708 

46.4060 

4.2904 

4.8999 

.03199 

40.0905 

4.3225 

4.7999 

.03879 

33.7500 

4.3062 

4.6998 

.04933 

27.1034 

4.3338 

Section 

Thickness  400.00 

5.0000 

.02512 

50.0301 

4.9000 

.03111 

41.2115 

4.4079 

4.8000 

.03846 

34.0337 

4.3302 

4.6999 

.04883 

27.3780 

4.3458 

4.5999 

.06819 

20.0320 

4.3531 

4.4999 

.11357 

12.2945 

4.1959 

4.4000 

.15149 

9.4265 

4. 3 COO 

.17634 

8.2863 

3.5422 

*.2000 

.20128 

7.4324 

3.3111 

4.0999 

.22937 

6.6814 

3.3545 

3.9999 

.25763 

5.8694 

3.4094 

3.89^9 

.32599 

4.9423 

3.4869 

3.7999 

.43139 

3.8330 

3.4694 

3.6999 

.61350 

2.7681 

3.4616 

Section 

Thickness  140.00 

3.9000 

.32410 

4.9709 

3.9114 

3.8000 

.43142 

3.8327 

3.4763 

3.7000 

.61704 

2.7521 

3.4545 

T 

57.3296 

51.2881 

45.1023 

38.7565 

32.2275 

25.4277 

17.9272 


62.8706 

57.0707 

50.9751 

44.9442 

38.4821 

32.1265 

25.2612 


39.2169 

32.2991 

25.4449 

17.8824 

11.2612 

8.0520 

7.2619 

6.4906 

5.6467 

4.6556 

3.5446 

2.4623 


8.4005 

3.5371 

2.4586 
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TABLE  VI.  (Centr'd) 


SECOND  HIGHER  RAYLEIGH  MODE 


Section  Thickness  400.00 


c 

k 

4.9000 

.05537 

4.8000 

.06549 

4.6999 

.08015 

4.5999 

.10807 

4.5000 

.25165 

4.3999 

.29310 

4.2999 

.32739 

4.1999 

.36497 

4.0999 

.40722 

3.9999 

.47463 

3.8999 

.59153 

3.7998 

.84393 

T 

U 

23.1581 

4.2866 

19.9891 

4.2627 

16.6789 

4.2938 

12.6396 

4.3836 

5.5484 

3.9105 

4.8721 

4.4642 

4.4632 

3.2661 

4.0990 

3.1949 

3.7634 

3.4253 

3.3096 

3.4644 

2.7237 

3.5669 

1.9593 

3.4754 

Section  Thickness  140.00 


4 . 6000 

.09402 

4 . 5000 

.25218 

4.3999 

.29289 

4.2999 

.32734 

4.1999 

.36532 

4.0999 

.40970 

3.9999 

.47503 

3.8999 

.58812 

3.7998 

.84289 

14.5273 

5.5369 

3.8832 

4.8756 

3.5369 

4.4640 

3.3107 

4.0951 

3.3051 

3.7406 

3.3130 

3.3'-69 

3.4600 

2  7395 

3.5842 

1.9617 

3.4812 
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T 

22.2932 

19.1184 

15.6739 

11.0562 

5.3337 

6.2623 

4.3787 

4.0190 

3.6417 

3.1683 

2.5201 

1.8086 


5.3339 

4.7581 

4.3744 

4.0019 

3.6409 

3.1671 

2.5214 

1.8076 


TABLE  VI.  (Cont'd) 


MODEL  314 


FUNDAMENTAL  RAYLEIGH  MODE 

Section  Thickness  810  km 

c  k  T  U  T 

4.5000  .00800  174.4980  3.6546  152.3970 

4.4000  .00910  156.9458  3.6739  133.7163 

4.3000  .01054  138.6728  3.7362  111.4856 

4.1999  .01273  117.4773  3.8025  84.3855 

4.0999  .01661  92.2917 

Section  Thickness  310 


k  T  U  T 


4.0000 

.02547 

61.6759 

3.4695 

39.4730 

3.9000 

.03856 

41.7795 

3.1635 

33.1033 

3.8000 

.04824 

34.2777 

2.9694 

29.2182 

3.7000 

.05646 

30.0793 

2.3125 

26.2169 

3.5999 

.06486 

26.9112 

2.7168 

23.4925 

3.4999 

07445 

24.1143 

2.7635 

20.7638 

3.3999 

.08602 

21.4831 

2.8034 

17.8086 

3.2999 

.10206 

18.6555 

2.8919 

13.9951 

3.1999 

.12901 

15.2204 

FIRST  HIGHER  RAYLEIGH  MODE 

Section 

Thickness  310  km 

c 

k 

T 

U 

T 

5.0000 

.02685 

46.8048 

4.4483 

36.2523 

4.9000 

.02886 

44.4386 

4.3098 

27.9880 

4.8000 

.03353 

39.0360 

4.2401 

21.5914 

4.6999 

.04475 

29.8741 

4.1392 

16.2429 

4.5999 

.05888 

23.1967 

3.8165 

13.1554 

4.4999 

. 08004 

17.4449 

3.5115 

11.5744 

4.3999 

.10431 

13.6899 

3.3171 

10.4853 

4.2999 

.12285 

11.8943 

3.2139 

9.5683 

4.1998 

.13938 

10.7339 

3.1330 

8.6651 

4.0993 

.15650 

9.7925 

3.1344 

7.7859 

3.9999 

.17704 

8.8727 

3.1097 

6.8751 

3.8999 

.20128 

8.0045 

3.1614 

5.9217 

3.7999 

.23306 

7.0948 

3.6999 

.27476 

6.1807 
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TABLE  VI-  (Cont'd) 


MODEL  3.1  j 


FUNDAMENTAL  RAYLEIGH  MODE 


Section  Thickness  810  km 


c 

k 

T 

U 

4.5000 

.00820 

170.2393 

4.4000 

.00930 

153.5283 

3.6992 

4.3000 

.01086 

134.5723 

3.7313 

4.1999 

.01334 

112.1333 

3.8018 

4.0999 

.01812 

84.5767 

3.8692 

Section 

Thickness  310  km 

c 

k 

T 

U 

4.0000 

.03053 

51.4482 

3.9000 

.04399 

36.6250 

3.4296 

3.8000 

.05404 

30.5995 

3.2226 

3.7000 

.06472 

26.2377 

3.0726 

3.5999 

.07614 

22.9218 

3.0533 

3.4999 

.09221 

19.4693 

3.0333 

3.3999 

.11718 

15.7708 

3.0416 

3.2999 

.16256 

11.7130 

3.0354 

3.1999 

.25847 

7.5969 

3.0323 

FIRST  HIGHER 

RAYLEIGH  MODE 

Section 

Thickness  310  km 

c 

k 

T 

U 

5.0000 

.02700 

46.5387 

4.9000 

.02898 

44.2546 

4.8000 

4.6999 
4.5999 
4.4999 
4.3999 

.2999 
.1998 
.0998 
.9999 
3.8999 
3.7999 

3.6999 
3.5998 
3.4998 
3.3998 


.03375 

.04593 

.06249 

.09907 

.12798 

.14814 

.16823 

.19220 

.22623 

.27637 

.36302 

.49743 

.72066 

1.18182 

3.20689 


38.7906 

29.1047 

21.8596 

14.0942 

11.1578 

9.8638 

8.8928 

7.9737 

6.9435 

5.8296 

4.5549 

3.4140 

2.4220 

1.5191 

.5763 


4.4580 

4.3170 

4.3459 

4.1502 

3.7883 

3.4936 

3.4019 

3.4011 

3.4963 

3.4527 

3.3824 

3.3589 

3.3175 


T 


148.7251 

129.1714 

105.2905 

74.8110 


T 


34.8049 

29.4172 

25.3270 

21.9835 

18.5069 

14.7044 

10.5790 

6.3232 


T 


35.9338 

27.2273 

19.5919 

13.0875 

10.7461 

9.6058 

8.6584 

7.7268 

5.4822 

/  r 

ZJi; 

3.1453 

2.1583 

1.2885 
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TABLE  VI.  (Cont'd) 


MODEL  310 


FUNDAMENTAL 

RAYLEIGH  MODE 

Section 

Thickness  400  km 

c 

k 

T 

U 

4.0000 

.02439 

64.4070 

3.9000 

.04255 

37.8649 

3.5912 

3.8000 

.05805 

28.4813 

3.2801 

3.7000 

.07023 

24.1816 

3.1144 

3.5999 

.08334 

20.9422 

2.9759 

3.4999 

.09981 

17.9868 

2 . 9384 

3.3999 

.12365 

14.9456 

2.9835 

3.2999 

.16535 

11.5152 

3.0190 

FIRST  HIGHER  RAYLEIGH  MODE 

Section  Thickness  140  km 

c 

k 

T 

U 

4,4000 

.12476 

11.4457 

4.3000 

.16172 

9.0352 

3.6088 

4.2000 

.18471 

8.0990 

3.3507 

4.0999 

.21052 

7.2797 

3.2917 

3.9999 

.24138 

6.5077 

3.4153 

3.8999 

.29117 

5.5333 

3.4070 

3.7999 

.36986 

4.4707 

3.4199 

3.6999 

.49765 

3.4125 

3.3786 

MODEL  320 

FUNDAMENTAL 

RAYLEIGH  MODE 

Section  Thickness  310  km 

C 

k 

T 

U 

4.0000 

.02676 

58.6913 

3.9000 

.03996 

40.3191 

3.4514 

3.8000 

.04957 

33.3568 

3.1758 

3.7000 

.05820 

29.1794 

2.9781 

3.5999 

.06716 

25.9897 

2.8662 

3.4999 

.07764 

23.1217 

2.7526 

3.3999 

.09056 

20.4061 

2.8481 

3.2999 

.11217 

16.9749 

2.8878 

3.1999 

.15391 

12.7579 

3.0158 

T 


34.7723 

27.2310 

23.2675 

20.2145 

17.2780 

14.1014 

10.4775 


T 


8.7490 

7.9020 

7.0921 

6.2579 

5.2734 

4.1854 

3.1475 


T 


38.1989 

32.1857 

28.3310 

25.2547 

22.4878 

19.5889 

16.0079 

10.8511 
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TABLE  VI.  (Cont'd) 


FIRST  HIGHER  RAYLEIGH  MODE 


Section  Thickness  310  km 


c 

k 

T 

U 

T 

5.0000 

.02694 

46.6523 

4.9000 

.02888 

44.4016 

4.8000 

.03358 

38.9835 

4.4617 

36.0753 

4.6999 

.04490 

29.7720 

4.2997 

27.9473 

4.5999 

.05957 

22.9299 

4.2645 

21.2086 

4.5000 

.08333 

16.7565 

4.4000 

. 11006 

12.9742 

3.8455 

12.4357 

4.3000 

.12959 

11.2758 

3.5636 

10.9456 

4.1999 

.14771 

10.1280 

3.3148 

9.8943 

4.0999 

.16573 

9.247C 

3.2854 

9.0110 

3.9999 

.18724 

8.3893 

3.1800 

8.1781 

3.8999 

.21377 

7.5367 

3.1997 

7.3062 

3.7999 

.24984 

6.6182 

3.2181 

6.3649 

3.6999 

-30270 

5.6103 

3.2358 

5.3289 
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from  6-10  sec.  Group  velocities  recorded  at  Oxford  indicate  a  slightly 
shorter  period  for  the  minimum  and  a  gradual  decrease  in  velocity  foi 
periods  shorter  than  12  sec.  However,  in  general  the  curves  are  quite 
similar  to  those  for  Atlanta  and  Blacksburg.  Waves  arriving  at  Atlanta 
from  Jalisco,  Mexico,  (parallel  waves)  have  both  fundamental  and  first 
higher  Rayleigh  mode  group  velocities  considerably  lower  than  do  the 
perpendicular  wave  trains  for  periods  less  than  about  20  sec.  The  rapidly 
"tailing-off"  of  group  velocities  below  15  sec  is  probably  due  to  a  thick 
sedimentary  sequence  (approximately  400  km  of  the  travel  path  for  these 
waves  lie  across  the  Gulf  Coastal  Plain). 

The  most  significant  difference  is  between  the  Blacksburg  and  Atlanta 
group  velocities  for  parallel  wave  trains.  A  minimum  of  2.75  km/sec  occurs 
at  a  period  of  22  sec  and  a  maximum  of  2.80  km/sec  occurs  at  a  period  of 
16  sec  for  the  Blacksburg  velocities.  Group  velocities  at  periods  greater 
than  30  sec  trend  toward  those  for  the  perpendicular  waves.  First  higher 
mode  Rayleigh  group  velocities  for  Blacksburg  are  lower  by  about  .1  km/ 
sec  than  those  for  Atlanta. 

Several  models  were  constructed  using  the  variation  of  phase  velocity 
with  layer  parameter  curves  for  the  Gutenberg-Birch  II  model.  Basic 
differences  in  the  models  are: 

1)  314,  315,  and  320  have  a  crystal  thickness  of  50  km,  while  310 
and  G-B  (Gutenberg-Birch  II)  have  a  crustal  thickness  of  40  km; 

2)  314  has  a  low  velocity  zone  centered  at  15  km  in  the  upper  crust. 
Except  for  this  low  velocity  zone,  320  is  identical  to  314; 

3)  G-B  and  310  have  a  low  velocity  zone  in  the  upper  mantle  beginning 
at  60  km  while  the  low  velocity  begins  at  70  km  for  314,  315,  and  320. 

4)  Velocities  in  the  first  10  km  of  G-B  are  slightly  higher  than  those 
in  the  other  models. 
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For  waves  traveling  perpendicular  to  the  Appalachians,  the  group 
velocity  data  agree  well  with  model  310  or  G-B  values  for  the  fundamental 
Rayleigh  mode  at  periods  less  than  30  sec.  At  longer  periods,  the  values 
fall  below  those  for  310  and  approach  those  for  314  and  315.  Fi’^t  higher 
Rayleigh  mode  values  observed  at  Blacksburg  lie  between  those  for  314  and 
315  for  periods  greater  than  10  sec.  Thus,  it  appears  that  waves  arriving 
at  Atlanta,  Oxford,  and  Blacksburg  from  the  Southern  Alaska  epicenter 
are  just  beginning  to  "feel"  the  Appalachian  structure.  The  length  of 
the  "Appalachian  path"  is  about  240  km  in  the  Appalachian  foreland  in 
each  case, assuming  that  the  Appalachian  structure  extends  beneath  the 
Mississippi  Embayment  (Oxford  station). 

Waves  arriving  at  Atlanta  from  Jalisco,  Mexico  yield  group  velocities 
considerably  below  those  of  G-B  for  periods  less  than  20  sec  (See  Fig.  11) . 
Fundamental  mode  group  velocities  at  Blacksburg  show  fair  agreement  with 
the  theoretical  curves  for  either  models  314  or  320  being  considerably  below 
curves  for  G-B,  315,  and  310  at  all  periods.  On  the  basis  of  fundamental 
mode  group  velocities  it  is  impossible  to  determine  which  of  models  314 
or  320  most  closely  approximate  the  crustal  structure.  However,  first 
higher  Rayleigh  mode  group  velocities  from  model  314  give  a  considerably 
better  fit  to  the  observed  data  than  do  those  from  v^oiel  320  in  the  period 
range  from  6  to  15  sec.  In  this  period  range,  the  group  velocity  curve 
for  314  lies  between  that  for  320  and  the  observed  group  velocities.  Thus, 
a  slight  velocity  reversal  in  the  crust  is  indicated  by  the  first  higher 
mode  group  velocities. 

An  alternative  to  the  low  velocity  zone  is  to  lower  the  velocities 
and/or  densities  in  the  first  10  km  of  the  crust.  However,  refraction  data 
give  a  compressional  velocity  of  5.88  km/sec  for  the  upper  crust  which  has 
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been  used  in  models  314  and  320.  The  shear  velocity  value  of  3.88  km/sec 
in  models  314  and  320  corresponds  to  a  Poisson  ratio  of  .25. 

Group  velocities  for  periods  greater  than  about  30  sec  indicate  a 
mantle  low  velocity  zone  beginning  at  about  70  km. 

Compressioral  velocity  crustal  and  upper  mantle  structure  to  the  west 
of  the  Southern  Appalachians  determined  from  the  travel  times  of  local 
earthquakes  and  for  the  Southern  Appalachian  structure  as  determined  from 
Rayleigh  wave  dispersion  are  shown  in  Fig.  12.  Gravity  and  Rayleigh  wave 
dispersion  data  indicate  a  total  crustal  thickness  of  about  50  km  beneath 
the  Southern  Appalachians.  Travel-time  and  dispersion  data  indicate  an 
upper  mantle  velocity  of  8.10  km/sec.  Dispersion  data  Indicate  a  slight 
low  velocity  zone  in  the  upper  crust  and  a  general  increase  of  velocity 
and  density  with  depth  below  this  zone.  With  the  exception  of  the  higher 
crustal  velocity  in  the  first  2  km,  the  Southern  Appalachian  s  ture 
approximated  by  model  314  is  similar  to  the  Northern  Alpine  structure  with 
a  50  km  crust  reported  by  Knopoff  et  al  [1966]. 
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7.0  CONCLUSIONS 


The  conclusions  which  have  been  drawn  from  this  study  are: 

1)  Focal  depths  of  local  earthquakes  in  the  Southeastern  United  States 
are  shallow  ranging  from  about  7  to  18  km. 

2)  Systematic  deviations  in  P-residuals  observed  at  Chapel  Hill, 

North  Carolina,  and  McMinnville,  Tennessee,  have  magnitudes  from  +  3  to 
-3  sec.  The  deviations  indicate  systematic  errors  in  the  Jeffreys-Bullen 
travel-times . 

3)  Travel-time  curves  constructed  from  local  earthquakes  and  refraction 
data  indicate  a  crustal  structure  for  the  Appalachian  foreland  of  h^  - 

33.0  km  (a  *  5.88  km/sec),  *  10.8  km  (a  ■  6.58  km/sec),  and  an  upper 
mantle  velocity  of  8.10  km/sec. 

4)  Rayleigh  wave  dispersion  data  Indicate  a  crustal  low  velocity  zone 
centered  at  about  15  km,  an  upper  mantle  low  velocity  zone  beginning  at 

70  km,  and  a  total  crustal  thickness  of  50  km  beneath  the  core  of  the 
Southern  Appalachians. 

5)  Gravity  data  indicate  a  total  crustal  thickness  beneath  the 
Southern  Appalachians  of  at  least  50  km  which  is  in  agreement  with  the  Rayleigh 
dispersion  data.  Gravity  data  indicate  a  crustal  thickness  of  about  45  km 

in  central  Kentucky  thickening  eastward  to  about  50  km  beneath  the  core 
of  the  Appalachian  and  thinning  to  23  km  beneath  the  North  Carolina 
continental  margin  at  about  the  2400  fathom  contour. 

6)  Higher  Rayleigh  modes  can  be  observed  using  digital  filtering 
techniques.  Care  must  be  taken  to  Isolate  Rayleigh  type  motion  and  to 
remove  interference  effects.  This  is  particularly  true  for  higher  modes 
than  the  first. 

7)  It  should  be  possible  to  use  higher  modes  than  the  first  to  delineate 
fine  detail  iu  the  crust. 
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APPENDIX  I 

EVALUATION  OF  THE  SECULAR  EQUATION  IN  COMPUTING  RAYLEIGH  DISPERSION 


In  the  Dunkin  Method,  the  Rayleigh  dispersion  equation  is  the  secular 
equation 
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Writing  out  (1-1)  explicitly  yields 
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where  the  ...  signifies  multiplication  of  the  form  g  p  ^ ...  g  *|”2 
Continuing  the  process 

n  i.  d  P-1  P-2  P-2  112  .  p-1 112  p-2  113 

Det  R1X  -  dxg  |J2  g  ...  +  axg  |13  g 


+  4  more  similar  terms  with  e,  as  a  factor 

+  5  more  expressions  with  a2,  a^,  a^,  a^,  and  a^ 
replacing  a^ 


where 


(1-4)  is  o£  the  form  (1-2) ;  the  sum  of  six  terms  which  are  each  products 
of  the  second  order  subdeterminants.  The  process  of  evaluating  the  second 
order  subdeterminants  of  a  given  layer,  multiplying  by  six  previously 
determined  coefficients,  s4 ,  and  summing  to  obtain  six  new  coefficients, 

,  is  repeated  until  the  entire  layered  sequence  has  been  traverseo  from 
bottom  to  top.  This  can  be  quickly  and  easily  adapted  to  computer  processing 
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as  is  indicated  in  the  following  considerations  of  the  real  and  imaginary 


structure  of  1 1“  and  g  P  ^ 


is  of  the  form  of  i  6  s  i  row  matrix 
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T..o  .real  czid  imaginary  structure  of  the  g|r^  is  of  the  form  of  a  6  x  6 
mati'-x 
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The  ab  indices  hive  been  taken  as  the  row  indices  and  the  cd  indices  the 


column  indices  with  1-12,  2-13,  3-14,  4-23,  5-24,  and  6-34; 
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Multiplication  of  tj^  by  the  gF  in  (1-2)  yields  a  matrix  of  the  form 
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Elements  of  this  matrix  correspond  to  the  terms  a.jb^,  ^^l*  etc"  *n 
(1-4).  Thus,  the  six  new  coefficients  A^  are  c ’Imply  the  su*"  uf  the  elements 
in  a  column  of  (1-7).  The  air  new  coefficients  can  then  be  considered  as  just 
another  6  /.  1  row  matrix  of  the  form  of  (!•  :  >,  ftinus  ones  have  been  inserted 
in  (I-f.)  to  account  for  the  multiplication  of  two  like-signed  imaginary  quantities. 


APPENDIX  II 


MODAL  SHAP h  COMPUTATIONS 


In  the  Dunkin  Method,  mode  shape  Is  computed  from  the  relation 


Rnm(z;a) 
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where  a  denotes  the  component  •»£  displacement.  R^  i.z;a)  is  the  a^ 
component  of  the  displacement-stress  vector,  at  the  depth  z,  normalized 
to  the  vertical  displacement  at  the  free  surface.  Note  that  r^.,^-  is  simply 
a  constant  multiplier  for  all  of  the  displacement-stress  components.  Thus, 
it  may  be  dropped  and  the  Rnm(z;a)  at  different  depths  divided  by  the  vertical 
displacement  rt  the  surface  to  normalize  the  displacements.  The  equation 
actually  used  for  the  determination  of  mode  shape  was 
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Writing  (II— 2)  explicitly 
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When  the  dots  denote  the  product 
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Now  the  g1^  which  have  explicit  subscripts  can  be  included  in  the  constant 
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factor  a^.  Collecting  ti  rms  with  common  factors  of  g^  yields 
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which  is  of  the  form  of  (II— 3) . 


Continuing  this  process  yields 
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The  Horizontal  component,  u  ,  of  displacement  corresponds  to  a*l  and  the 
vertical  component,  w  ,  to  a*2,so  that 
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Writing  (II— 7)  explicitly 


u  ■  a 


n 

12 

n-1 

ji. 

n  112  n-1 

I13 

n 

I12 

n 

a2 

g 

12 

g 

lef 

. . .  + 

a2 

8  |l3  8 

|ef  *•* 

+  a2 

g 

I14 

g 

n 

12 

n-1 

|23 

n  |12  n-1 

24 

n 

112 

n- 

a2 

g 

23 

g 

lef 

...  + 

a2 

8  |24  8 

ef  *•* 

+  a2 

g 

I34 

g 

n 

13 

n-1 

112 

a3 

g 

12 

g 

lef 

...  + 

o  more  terms  in  a^ 

n  | 

14 

n-1 

112 

in  a^ 

a4 

g 

12 

g 

lef 

...  + 

5 

more  terms 

u  is  expressed  in  terms  of  the  sum  of  18  terms  of  the  fora  a  gn  gn  ^ 

m  |lci 

Collecting  terms  with  common  gn  *  gives  u  as  the  sum  of  six  terms  of 
the  form  am  gn-1  . . . 


,  n-1 112  .  ,  n-1  Il3  .  ,  n-1 114 

a  1  8  ef  ”  ’  +  a  2  8  |ef  *  ”  +  a  3  8  |ef  * '  *  + 


,  n-1 123  .  ,  n-1  124  .  ,  n-1  134 

1  4  8  |ef  a  5  8  |ef  •**  +  a  6  8  |ef  *  *  *  » 


(II-7) 


(II— 8) 


(11-9) 


(11-10) 


II-3 


where 


“1 


t 


«n  (12  +  ‘3  Is 


etc. 


The  process  of  evaluating  the  displacements  u  and  w  now  proceeds  exactly 

as  for  the  evaluation  of  Det  described  in  Appendix  I. 

Let  us  consider  the  real  and  imaginary  forms  of  the  factors  appearing 

in  (II-4),  (II-6) ,  and  (II-1Q).  From  the  definition  of  ^  given  by 

Dunkin,  t,  is  of  the  form 
Ir 


t^r  ■  [I  R  R  I]  ,  I  ■  pyre  imaginary  quantity 

R  ■  pure  real  quantity 


(11-11) 


Tiie  laver  matrices  G  are  of  the  form 

n 


G 

n 


R  I  I  R 
I  R  R  I 
I  R  R  I 

mu 


(H-12) 


Thus  t.  gP_1 
lr  ®rs 


[I  R  R  I] 


(11-13) 


Each  successive  multiplication  by  a  layer  matrix  G  yields  a  row  matrix 

n 

of  the  form  [I  R  R  I]  until  the  first  multiplication  by  a  second  order 
sub determinant  is  encountered.  In  order  to  account  for  the  multiplication 
of  two  imaginaries  in  the  evaluation  of  column  2  and  3  in  (11-13)  minus 
ones  are  inserted  in  (11-12)  to  give 


II-4 


*  v  * 


For  the  u  component 


1 


[R  R  I] 


R  R  I  I  R  R 
R  R  I  I  R  R 
I  -I  R  R  -I  -I 


[R  R  I  I  R  R] 


and  for  the  w  component 


[I  R  I] 


R  R  I  -I  R  R 
II  R  R  I  I 
R  R  -I  -I  R  R 


[I  I  R  R  I  I] 


(11-17) 


r 


(11-18) 


Minus  ones  have  been  inserted  to  account  for  multiplication  of  two  imaginaries. 


To  continue  the  process  to  the  surface,  the  form  of  the  matrices  composed  of 
the  second  order  subdeterminant  elements  is 


II-5 


R  R  I  I  R  R 

R  R  I  I  R  R 

-I  -I  R  R  -I  -I 

-I  -I  R  R  -I  -I 

R  R  I  I  R  R 

R  R  I  I  R  R 


(11-19) 


for  the  u  component  and 


R  R  -I  -I  R  R 
R  R  -I  -I  R  R 
II  R  R  I  I 
II  R  R  I  I 
R  R  -I  -IX  R 
X  R  -I  -I  R  R 


(11-20) 


and  for  the  w  component.  Minus  ones  have  been  inserted  in  the  appropriate 
elements  to  account  for  the  multiplication  of  two  pure  imaginary  quantities. 
Since  the  lower  subscripts  in  (Il-l)for  the  surface  layer  matrix  are  21,  the 
layer  matrix  of  the  surface  layer  is  a  6  x  1  column  matrix.  The  final  multi¬ 
plications  are  then  of  the  forms : 


u  component 


[R  R  I  I  R  R]  * 


[  R  R  R  R  R  R] 


H-6 


EXPLICIT  FORMS  OF  LAYER  MATRIX  COMPONENTS  g 


SUMMARY  OF  EXPLICIT  FORMS  OF  SECOND  ORDER  SUBDETERMINANTS  g  |j^ 


Sl2  *  g‘  34  "  -  2.Y(Y  -  1)  +  (2y 2  -  r  +  DCP  CQ  -  +  Y^rgjSQ  SP 


112  24  ,  ,2,v-l  ,SO  CP  - , 

|13  “  8  34  "  4  (pX  k)  I  r -  +  ra  sp  CQ  ] 

B 


lu  “  8  23  a  8  |34  "  8  [34  “  Up*2*)'1  |(2y-l)  (1-CP  CQ)  +  +  yr^SP  SQ 


2  13 


-  (pX2k)-1  [r  CP  SQ  +  —  SP  CQ  ] 


8  £4  -  -  (P*2k)~2  [2 (1-CP  CQ)  +  (•— -+  r  rj  SQ  SP  ) 

H  rar0  a  15 

8  |l2  “  8  (24  “  "  (Px2k)  [Y2rg  CP  SQ  +  SP  CQ  ] 

8  Is'*  £4-^® 


II-9 


I1*  ••  |23  |34  134  .  ,2,  f  ,  1W. 

8  ( 12  "  s  Ji2  "  8  1 14  *  8  U3  "  i(pkX  M  Y(Y-1)  (*-l> 


+  [■^L=^"-  +  r  r.y3]  SO  Spi 

Vs  a  8  > 


|14 

123 

g 

|l3  "  8 

13  ’  8 

14 

23  ,  t 

g 

14  *  8 

23  “  1  4 

14 

23 

g 

23  *  8 

14  “  8 

124 

34  A  , 

8 

12  “  8  | 

13  *  + 

24  r 

g 

13  *  +  r 

-  SP  sc 

B 

24 

14  *  123 


8  cp  SQ  +  r 

S' 


ex' 


|3  -  1  +  2y(y-1)  (1-CQ  CP)  +  (  +  y2 


rar8 


14 

14 


-  1 


It  -  +  (pX2k)  CP  SQ  +  rY2  SP  CQ 


1*34  2  2  2  2  —  — 
g  P  -  (pAZk)  i2y/(y-1)  vCP  CQ  -i) 


r(y-l)  .  4 

[  _  _  +  Y  I 


Vs 


11-10 


(1-CQ  CP) 


CQ  SP  ) 


Vs)  SPSQ 


] 


r J  SP  SO 
a  p 


11-11 


SP  ■  siuh  P 
CP  ■  cosh  P 


SQ  ■  sinh  Q 
CP  «  cosh  Q 


c  <  a 


c  <  a 


The  expressions  are  written  for  the  hyperbolic  functions  (c<ct,  c<3) . 
For  Ooi,  c>8  ,  the  hyperbolic  functions  transform  according  to 


C"8h  P 

- >  cos  P 

c>a 

cosh  Q 

- »  cos  Q 

OB 

sinh  P 

- >  -1  sin  P 

c>a 

sinh  Q 

- >  -i  sin  Q 

c>B 

11-12 


APPENDIX  III 


FORTRAN  LISTINGS  OF  COMPUTER  PROGRAMS 
FLAT RAY ,  STRESS,  INTEGRAL,  TRAVEL,  VARGRAV,  AND  PRESID 


INPUT 

FORMAT 

14 

14 

14 

F10.4 

FlO.4 

8F10.6 

F10.6 

F10.6 

F10.6 

F10.6 

F10.6 

FlO.f 

F10.6 

F10.6 

?0I4 


DATA  AND  FORMAT  FOR  FLATRAY 

DATA  DESCRIPTION 

Model  Identification  no. 

No.  of  layers  <41 

No.  of  modes  to  be  computed  <11 

Precision  of  root  values  of  k 

Cutoff  k  value  to  automatically  stop 

computation 

Layer  parameters  stored  In  sequence 

a*  8,  ,3,  d;  two  layers  to  a  card 

Starting  c  value 

Minimum  c  value 

Decrement  of  c 

Starting  k  value 

k  Increment  for  individual  modes 
k  increment  to  find  starting  points  of 
different  modes 
c  perturbation 
k  perturbation 

Mode  no.  to  be  found  stored  sequentially 


III-l 


AlPha  mo  KBfcMAUO)  ,  HhO  <  4  0  ) .  ThI  CK  *.  40  >  .ROOT  St  15,40  )  .MODE 
It lUi .PERIOD < 15.40  > ,dl 6 ),£<*> ,G< 4.6) , At  6>,UDISP< 40) ,*UlSPt40 ) 

1.M44I 

CQHNU*  UOISW.»»OlSr  ,h 

«=AU  15Q.H0UEL.NLAUR, RHODE. X*HJN  , 4KHAX 
150  ► 3«HA  T  ( 3M.2F10  .4) 

He AQ  200. < ALPHA  I I  )  ,  HgT  T  A (  I  ),RhU( I  ) .  TH J CK t  I ) , I « 1 . *1 A ffc R ) 

200  ► }*HAT  tpFl 0 .6) 

PRU*  900-HUDtL 

900  T  DHiAT  t  01 H* LAYER  PaHaHETC-RS  fO«  HOD-‘L  RTI  14  //  6H  LAYER,  5X,*H  AL 
lP4A.9A.6H  BSTTa.7A.4H  HH0.5*»14H  T H  I  CK NESS ( K H )  ,  3 X  ,  6H  H<KH>  /> 

I T* \lt *0  . 

UD  9U2  I«1,NLAYE« 

»  ThicpuTThick  ♦Thick  1 1 ) 

P-UN»  902,  1 .alpha t I ),B6TTAt  P.HHOtl ),ThICK<I),TThICK 
9y2  ►  3HMAH1A.  t3.OX.FlO.4.4X.f0.4.4X.FB.4,5A,rS.?,5X.FB.2) 

I  )  I  AL  «  TThICK 

«iAO  200.  CHAX.cHIN,DfcLC»AK,D6LK,DbwKl .PTRBC.PTRB  K 
h • a 0  J5  0 ,  < NODE < I ) . I >1 .RHODE ) 

150  ► 3*4aT  <  20 l 4 ) 
laOM*0 
NiHHUO  «  1 
I'MAXl  s  CH  A  X 
*400  i  10 
t 4«HU0b< 1 ) 

151-41  >10 
U>  nlayeh-i 
l  4«1 

A4S4U0  lit  I  H  ,  INMAK 
I  3  I  4b  >1 
l  s  CH A  X 
Jj^i«l 
500  JOJS^s  1 
jjHa«l 
13P1*1 

UJI|R 
i  I  vi9«  1. 

1 34 AC  >  10 

i 9UX2S1 

10 0  I  =  t  A 

UHC 3A*XX*C 

HAu^H  >  SiRHABSt  (C/ALPHA  (  I  X  ♦  1  .*  >**2-1.  > » 

"dcMSUHTi A8S( CC/bETTAf I  X ♦ 1 » ) • • 2- 1 . ) ) 

A.«2. •AK«*2-0HfcGA»«2/tBETTA( |X*1>**2> 
h •>  aL* *  bET  T  a(  |X*1  )**2/<  2.*0HEGA**2) 

AH0*»EITAC I  X ♦ 1 | • • 2*HH0 t 1 X ♦ 1 > 

♦HI  AC*AHU*<-4.«AH#»4**AlP4*HBEWAl**2> 

d( i  )  >Tl2«dErAC»RALPH«XK*<2.*Xx**2-A. ) 

t*<  3>«T13s0EFAC*XK*(-2.*XK**2*RAl*H*48ET*Al> 

H(4)*Tl4sTlO 

d<*»)*Tl5s-^dET»Tl2/HALPH 

d( 0)>Tl6>-BeF AC»XR**2»(-HALPH«HdET*l . )/AHu 
>1  B  3  T  T  A  <  I  )  »*2/C**2 

U| v$«i . 

x-iu^sR.iUt  I  )*C»»2**K 

«.lPh>5UHT( ABSt <C/ALPHA( I ) ) 
h^cT  tSQRT  t  ABSt  (C/BETTAt  I  > )  •  •  2-  1 .  )  ) 
h.x<MaIP-iMh|CM  |  ) 

U»A<*R0ET*TmICK  t I ) 

I'C-ALPHAt  I  )  11,2,2 

2  ♦ ALA*1 . 

CasCUStP) 

5a*5IM(P> 
liD  TO  3 

l  y  ala  «-l , 

L'3*(txP(p)#E*Pt-P)  )  /  2. 
ia*<tAP<P)-CAPC-P))/2. 

3  l- tC-dETTA< I) >5.4,4 

4  y  al;»si  . 

LJsCu'ilO) 
i 3  *S ! N ( u ) 
u3  TO  6 

?  y  ALd*- 1 . 

C3> ttXPtOl*EAPI-U) )/2. 
ijsi tAPIO  >-EXP< -0) 1/2, 

UJ  49«CP»C0 

o  l# »  i  .  1  )I  IGAH-1 .  U«2/tNAlPH*RB6T  )  «GAH**2*HALPH*RBET*F  AvA*FACB 

l.1 1 , 1  )*-G< 1 . 1 )»SQ*SP-2, *CAH«t GAH-1. ) ♦ ( 2, *GAH**2-2, *GAH«1 . )»CP*CO 
b(  1.2)ilS3»C;P/HflET*HALPH*SP*C0«FACA  »/«HOC 

t*(i,i)iU3A--l.  )/tHALPH#RBET)»UAH«KALPH»RBET*FACA«f  ACB)*SP*S0»«2.* 
’.OArt-1,  )*t  l.-CP»COl 
bt 1 , 0) >Gt 1 . 0) /HHOC 

b(  1  ,«>■(,  (  1.0) 

U<  1 ,9)*»<*8ET»CP*SG*f  ACB*SP*r  J/RalP-!>/RhOC 

o<  1 . 6)  «-  <  2.  •  <  1 . -CP»CU  >  ♦  t  1  ./<  -(ALPH»Hd6T  )*RaLPh«RBET»FACB»FaCa  )»su*s 
lH)/(>  -MJC**2) 

M2.1  )*.RhOC*<GaH**2*R8ET«FaCB*CP*»J*<GaH-i .  >**2*SP*CG/RaIPh > 

Ul 2.2)«C0»CP 

b( 2 . 0  )  «GAH«H8tT •FaCB*CP*S0* t GAH-1 . ) *SP* CO/ RALPH 
b(  2U)*bl  2.3) 

bt  2.5)sHBET*S0«SP*FaCB/RaLPh 
bt  2  .M*b<  1 .5) 

b( 0. 1 ) « ( ( 3AH-1 . )##i/ <KaIPH.RBE r ) ♦Ra.Ph»RBE T *GAH* • o*F AC8*FaCa )*SU 

1*»p 

b<3.1)«-  4H'*;#(GAH«tGAH-l. )*<2.*GAH-1.)«(1.-C0*CP)  tG(i,l)) 
b( i.2)>  <  tuA4-l.  )*CP*SO/HBET ♦RALPH* jAH*FACA*CO*SP) 

b!  1,0)11. ♦2.*b*H#t GAH-1, )#t l.*LP*CO)*t (GAH-1, ) * . 2/ t RALPH*RBE T ) *GAH 
1*#2*haLPH*H«ET*FaCA*FACB)*SP*SU 
b ( 0  ,  a  )  *Gt J.31-1. 
b( 0.9)S-Gt 2.3) 
bt J.6)s-g<1.3) 
bt A.  i)sG(3, 1 ) 
b( 4 , 2)*6t 0.2) 
b(  4 , 0 ) >b t  3 . 0  )  *1 . 
b(4.4)«u<3.0) 


T II-2 


li  (4 
UM 
4.45 
b]5 
U<5 
(.15 
v»i  5 
0]5 
U(0 
1<J 


U<6, 
G46. 
Ul  6,, 
U(6, 
uj  6,* 
i»I  6,  i 
I  r  i  1 1 

62 

1  r  4  1 

101 

03  3! 

SJMxl 

03  41 
SJN« 

40 

C3NT 
A|  J) 

30 

C0*f 
03  8 
tti  j) 

a 

COHT 
1*1- 
G3  TI 

402 

F  *6| 
1*346 
*  tF/l 

•  • 

i  •  »  1  ■ 
CJNPi 

603 

K  U 1 

60< 

1  MOXi 

Ml* 
ii«r 
*4*- 
03  Ti 

605 

A  •  * 
A?l* 

I  r  4  Al 

6C  ' 

U  4  XI 

606 

1  f  4  Al 

60V 

1*1X1 

613 

4*3  Ti 

703 

JJH6' 
*42x 
*4*  1 
1  0HAI 
i* 3  T> 

706 

I0HAI 

►  l*F( 

Ml* 

704 

1  0HAI 
I?  4  Xi 

614 

*4*X' 

*42*; 
U3  Ti 

>13 

>  2«r 
M<* 

M s  Al 
4*3  Ti 

606 

Kill 

70/ 

1  *  (  V 

61/ 

I  ••  4  Al 

623 

A<  S 

urn 

Kim 

1014 

A  4  *  Ai 

4,-J  Ti 

61V 

*41* 
►  l*r 
*4*  Ai 

1  r  4  X' 

306 

|F4  1 

3000 

KUi 

705 

*  2'F 
*42* 

30J 

►  1«F 
x  «  * 

"30t 
^“4  I1 

1  F  4  | 

Vjj 

1  4  *  |  1 

i  •*  4  r 

616 

►  1  *  F 
A  4  1  * 
A  4  *  Ai 
u3  Ti 

61a 

A4  **' 

UrL< 
A  4  H  J  i 
1  >  1  4 

5>«-Q4?.3) 

o)*«cu.j; 

1  )•  MGAH-1.  )«*2«CH»SQ/HB6  l •HAL* 4*GAM*GAH«5P*CG*F ACA  )*HhOC 
d)  *HA4.P4«5P*SU«F  ACA/H8ET 

3) *-  G43.2) 

4) i-GU.2> 

5) »U(2.2> 

6) *Gtl.2> 

1)*U3A*1-1 , 


)«i-4/«HAlPH»»<8fel  >*GA1'  •4*«ALPM»H«6T*f  ACA*f  ACH»*SP*S 

) 


1) «MH0C«*2»<2.*GAN#GAH*| 8AH-1. ) ••*• 4  CP* CO- 1 .  ) -0(0.1) 

2) *G(5.1 ) 

4) *-Q(3.1> 

5) «642, 1> 

0 ) • G I  1.1) 

UISP-D62. 62.8012 
-1)402.402.101 

U  Jb 1.6 

U. 

0  Ml. 6 

sOH*8(K >*Gl K.  j) 
iHUfc 
•SUN 
iNUfc 
J*1.6 
*A4 J)/Dl V$ 

!HUfc 

1 

U  51 

1)*G<1.1 >*B42)*G(2.1 )*B| 3>*G43*l)»B44)*G<4.1)»e45>*G45.1 MB<6) 

.1) 

U 1  VS 

aTKT-i)  603.301.603 

UT6  FIRST  HOOT  OF  EACH  HODt 

*0*2-1)605.604,695 

t  a  1 0 

A  4 

*MPEU*1 
O  100 
xgSlF  ) 
xgbtFl) 

►  -AH)  606.607.607 
HSIP-CD-AF)  608. 60V. 60V 
05IF1-F )-AFl)  608.609.609 
«s4  *K1-*K  )-XKMH)  308. 308. 610 
u  4 70i. 7Q4 ) , JUHP 

*2 

A* 

Xd5( I XK»F1-XK1*F2)/«F1-F2) ) 

1  a  1 

o  i ioo,60 ) .juhpi 

LtQ 

2 

*42 

CiO 

"-XHD  614.613*613 
*«XK«IN 

X< 

V  4)50.60) .JUHPl 
*< 

*-XAH|N 

u  I  100.60  ) .JUHPI 
»4AC* 1 ) 306 • 706. 707 
"00-1)617.618.617 
►-*Fl>  619.61V. 620 
XK-2.-0Ei.41 
*-l . 

P-0. >1014.1014. 

K*Ofc4.<l 


.100 


u  5U0 
M 

P  ♦ S | GN»  Ofe  L  K 1 
"-0. >1014.1014.100 
sTnT-iMOOO, 508. 3000 
nHAC-10)  3U0, >08,706 
1 

*41 


'  ) 


*851 ( XK?bF1-X*bF2)/IFi-, 

6|  I  M  .  I  N  )  •  X  K 

U0UX.INI*  42,«3.1416)/1  XX*C) 
a THT - 1 )  950.901.950 
*♦1 

“i-HOOEINMUUt )  >  615.615.616 

*4 

**0fcl<l 

u  ’00 

*  »  2  .  *aK41* 
lx  N00T8UH-l.il 
•'*  .01*bfLK| 

1*10 
*400  ti 

b*j  Tu  5U0 

J 3***1 

1 400*1 

14  *NOOE 1 1  NO 0  ) 


!II-3 


) 0  >  s^«10 

J  JHP 1  *  2 

M*HUDTb(lN  ,1) 

).*1 

*4rt|rta  .01*»* 

| jh^«l 

)  M*  1 

dQl)  J  I  «L *  I  X 

*te*2.»XK»*K-(XA*C)»»2/(BETT A( | 4*1) ) »»2 
4b^d*AL*AL 

*AL*«iSUHTUBS(  (C/AIPMAUX*1)  )»»2-l,  )  > 

AMU«bETT*( lx*l)*»2*HHU( IX*1) 

B(1)*BETTA(  1X*1)»»2/(*K»C»»2) 

0(2) *-AL»b€TTA( l X*1 >»»2/(2, «***»3»C»C»RalPM) 
b(4)«-BfeTTAt  1X*1)»»2/U,»ANU«  «  *K*C  *  ••2> 

B( 4 ) »  BE  T  T  A ( !X*1)**2/(2.*AHU»HALPH»(C»XK)»»2) 

6001  0010*6010  *8009 

BOO*  BcliiU«afcTTA<  1  )»»2 

*1U«  8£T  TA( 1 ) • *2*NHU( | ) 

H  4LPH  a  S(jR  T  (  A8S  (  (C/ALPSA(  1  >  )  ••2*1 .  )  ) 

H8tr«SQHT(ABS( (C/B£TTa< 1)1 ••2-1.  ) ) 

PaM*HALP4*lHlC«(  1 ) 

U«MMBfcT»THlCM  1  ) 

) *(C-ALPha( 1 1 )  8002 * 6003*6003 
6005  ►  XC A* 1  * 

C»*CUS(P) 

HiP  ) 

UD  TU  8030 
0002  ► ACAa-l. 

i;»«lfcXIHP)*£*P(-r )  1/2. 
i»«(tAP{p)-fcXP(-P) i/2, 

6030  U(C-dETTUl))  6004*0005.8005 
0005  ► aC6*i • 

Cj*CuS(U) 

* js$)N(Q) 
iO  TU  8006 
e004  .•  AC8«-1. 

Cj*(txP(Q)*EXP(-Q) )/2. 

( t  Q) -t*p( .Q)  )/2 . 

80 Oo  H 1 U  «  BfeTTM 1 >»»2/(XK»C»C> 

X.«2* •XA*AK-< IK»C)*»2/BETTA< 1 )»»2 
«( 1.1) *2. •6cTbU«CP/ <C*C)-AL*BHi*CU/xK 

u<l.  <*)*-(  XL  •B1G»SP/(RXlPM"XA  )*2.»Bl3*RBET«S0*FACB*XK) 

G(  l.dM-BlGMCP-CUWAHL 

U< it •)•» U*»P/( AHO»  MALPH)*BlG*SU»MyfT»r  ACB/AMU 
m<2. i  )*-2, •01G»RALPM»SP»rACA»*K-AL#8lG»SO/( IK*«BEI ) 

0(2,  i  )i-Xl«01G»CP/XM2.»B lG»Ca»XK 

t*(2.  Jl*H*uPH»BlG»SP»rACA/A«U^d)3»Sli/(RBET»AHU) 

o<  2*  4 ) a-G( 1 » 3 ) 

0<5.1)*2.»A1U»AL«BIG»(CP-CQ) 

G<3.2|a-A*1U.Al»AL«BlG»SP/<RxL»M*Xox*)-4,»AHU*RBET*SU*fACB»Bltt»XKa 

1*4 

b(  3  *  3 ) >G( 2*  2 ) 

G<3, 4)a  -3(1*2) 

G(4,l)**4,#*MU»*K»*K»RALPM*BlG*S*,»rACA-AMU**L*AL#BtG#SO/ 
1Ua«a<«RBET  ) 

G<4,2)t-G(3*l) 
ti(  4 , 3 )  a  «G( 2 »  1 ) 

U(  4,  •  )  a b( 1. 1) 

U3  Oti 0 7  J3l»4 
6J«aU , 

U3  0006  K *  1 • 4 

6jH«6JH*g(A )»G(R. J) 

0008  C3*f  1  NUfc 
X( J)a$U* 

800/  (O-aTlNl'k 

UD  8011  J«l*4 
fe(J)«  B(J) 

H( J) «0( J)«A( J) 

6011  l’J*TiMU£ 

I" *  lL-1100, 63.04 
03  I » I  - 1 

GO  To  0001 
0010  t<5)M<l> 
t<6)  M<2> 
fc(  /»*0(5) 
t (0)«a<4) 

)’(  1-1)81.61*02 
61  'OlOilHlCXt 1 ) 

HlCMI  )>0. 

!.>10 

!.*•*» 

G3  To  0009 

04  )?(  ltl-1)  65*05.06 

05  ■  4)C*(  I  )*H0i.D 

ui*$»i. 

GD  TU  51 

8012  I-  <  IU-DVO.07.00 

0/  U3l3H<2  )«-*(2)»G<l,l)-*(3)»G(2,l)-A(4)»0(3.1) 

■  0  1  S* ( 2  )«M1)»G(1,1MA(3)»G<4.1)-A(4)*G(5,1) 

uo)S»*(2)*  UD  1  SP(  2  )  /D 1  90 

"U 1 < 2 ) *  HDISP(2>/0196 

>,l«10 

H(1)«E( 1) 

B( 2 1 *E ( 2) 

0 <  3 ) >E ( 3 ) 

B(4)«E(4) 
go  ;g  6009 
00  l<4lC0(  1  )*0. 

U| 9S«1. 

G3  TU  51 

66  uuiSHi i  )a-A(2)*G( 1 . 1 ) -A < 3 ) *G ( 2, 1 ) -A ( 4 ) *R( 3 » 1 > 

■  QIS*M1  )*A(1)»G(  1. 1)*A(3)*G(  4, 1)*A(4)»G(  5,1) 

U0ISP(1)«  UOl SP ( 1 ) / D I  VS 

«D)S*M1)«  MU1SP<1)/01VS 


III-* 


lpisr*l 
UICM  I  )  *hOlO 
La  U*J 

►'JNCH  2701 ,ND DEI «  IM,L»T0TAL.C.PE81Dt)(  JM.  IN) 

►'JHCH  1032. (UDISP(U) »6DISP( Ui . Jal.D 
UIVlMOUP(l) 

UD  1001  J* 1 • I  * ♦ 1 
WO«SH( J)aJOISP< J)/DIV 
"OiSPC j)«KDISPt J)/Div 
1001  CONTINUE 

|3UN*|UUN*2 

PJNC"  2701,«00tL, IK. L» TOTAL, C.HBHIUjt |K. IN ) 

PjtfCM  1032. <UDISP( J) »NDISP( J), J«l.l) 

2701  ► DNK*T(3I4,3f9.3) 

1032  *  DKNAT (6E13.6) 

UO  TU  (8025. 901. 8026), IJKP 
82  >Olo*TmICM  I  ) 

I  N)  CM  I )  *  I) , 

i.*0 

UD  TU  800« 

89  UICM  I  )  **<010 
I L  I  *0 

Hal 

U|VS*1. 

UD  TU  81 

90  UD  Tu( 92,91. 93,96)  .  M 

92  UD  6014  K  « 1 , 6 
3  JH  a  0 . 

UD  8015  J«2.4 
3jM«*jM*B( J) *G ( J-l .* ) 

60)  >  CONTINUE 
A  (  K  ) a  SUM 
6014  CONTINUE 

UD  8031  J 1 1 , 6 
d( J) «A( J)/DI VS 
6031  CONTINUE 

)  a  |  -1 

M*2 

UD  TU  51 

91  |M  1*1 )  94.94,95 

93  UD  3016  0*1*6 
3J«*0. 

UD  6u 1 7  4*1,6 
SJN«»UK*B(K )»G(4 , J) 

601 7  CONTINUE 
A ( J ) a  SUM 

80 1  o  CONTINUE 

UD  6018  J« 1  *  6 
8( J>*A( JJ/OIVS 

6018  CONTINUE 
I  a  I  -1 

UD  TU  51 

94  UQ) Sr  It *1  )a-A(l)*G(l,l)*A(2)»G(2.1)-A(3)*C(3,l )-a(4)*C(4,1) 

l-A(5)aG*5.1 )*A ( 6 ) *G ( 6, 1 ) 

U0UiP(L«l)a  U0ISP(L*1>/DIVS 
b(l)«  M ( 1 ) 

8(2)*  M(2) 

8( 3) *M' 3 ) 

8(41*2(4) 

•  at 

MaO 

U I  V  5a i , 

UD  TU  51 

91  A ( 1 )  ■ -B ( 1 ) «G ( 1 , 1 )-B(3)*G(4,l) ♦8(4)»J(5,1) 
A(2)«-Bxl)*G(1.2)-B(3)*G(4,2)*8(4)#i(5r2) 

A(4)*3d )*G» A,3)*B(3)»GC4,3)-8(4)*G(5,3) 

A  (  4  )  a  t*  {  l  )aG(1.4)*Bl3)»GC4,4)«B(4)aU(5,4) 

A(5)»-B(l )*G(1,5)*B(3)»G(4,5) *B(4)»i(5,5) 
A(oJ«»B(1)*C(1.6)»B(3)*G(4.6)*H(4)*j(5.6) 

UD  6019  J* 1 . 6 
8  (  J  >  ■  A  (  J ) 

8019  CDN  ; NUE 
Ms4 

la  |-1 

UD  Tu  51 

9b  |M|-l)  97,97.96 
94  U( l,3)a-G( 1 .3) 

U<l,4)a-G( 1.4) 

U ( 2 , 3 ) a -G (2.3) 

U( 2 »  4 ) *-£( 2. 4 ) 

U( 3, 1 ) *-G( 3.1) 

U( 3 . 2 ) a-G( 3. 2) 

U(3.3)a-G(3.5) 

U(3.0)a-G(3.6) 

U ( 4 , 1  )  a-Q (4,1) 

U(4,2)a-G(4.2) 

UC4, 3)8-G(4, 5) 

U(4.6)a-G(4,6) 

Ul 9»3)s-Q( 5,3) 

U( 3  »  4 ) a-G( 5 • 4  ) 

U(  6 . 3 ) s-G ( 6 » 3 ) 

U(6,4)a-G(6.4) 

UD  6020  J  *  1 , 6 
3 JM 10. 

UD  8021  4*1,6 

tf(<)*G(4,j) 

6021  CD  if  I NUE 
A< j) * SUM 

8020  CDNT )NUE 

UD  6022  0*1,6 
6 ( J  > • A ( J ) /O I  VS 
8022  CONTINUE 
I  a  ) -1 
UD  TU  51 

9/  »0 1  SnL  •  1  )«-Atl)*GU.l)-A(2l»G(2.1)4A<3)*G(3.1)*AI4)*G(4,l) 


III-*, 


1-«45)*G<5, 1  >-A!6>*G!6.1> 

«OI»Nll*l>>  «DISP4L*1>/DIVS 
d(l)*E(5> 

UWXtift) 

dl3)sE<7> 

d(4)*fc!6> 

ltL 

L*L-1 

I  .** 

G3  *U  8001 

802?  Ma-«UDI*<  IN.  l)*D6t« 

I  W 

lM«MODt< 1HOU) 

Li  CHaX-DEIC 
2000  ►'■UN!  1030.1*.  TOTAL 

1030  ►  3**1*  T  < 6H-NQDE  13  /19H  SECTION  THICKNESS  F8.2  //11X.  15H  PHASE  V*L 
1  Or  I T » , 4*  # 2H  K,  hX.JOH  PERIOD!*) .4*. 15H  GROUP  VELOCITY, ex. 10m  PERIO 
1U<S) ) 

*MI4»  1031.CNAX.RDDTS! IN,l>.PeH|DD<  l*.l> 

1031  f 3*<*T! 121. r 10,4. 3X.F 10. 5,4*,F 10.4  > 

U3  Tu  5QQ 

301  irifNDXl-l)  502,501.502 

501  INUX2  ■  10 
*<1**4 
U»* 

*<  **4  ♦  OELK 
03  Tu  100 

502  A»  •  AbS ! F  ) 

A* l*A3S»!Fl> 

lflA>-Afl  t  504.503.503 

503  IMAdSiri-r)  -AD  507.506.506 

504  IDAdSIFl-F  >-Af  1  >507,506,506 
5Q6  I D AbS! X*l- *4) -XKK|N >508,508,610 
50/  IM|t»RAC-l>  508.706.513 

513  I*IA>-AF1>  510.510.511 
511  IDJUP-l)  2004,2004,510 

2004  IMJuPl-l*  2006.2006,2005 
200?  *<!**< 

►  l  *F 

u 3  lu  2006 
200?  »H*M«-1. 

*<•*4-2. •  DEL* 

1 3**  1*10 

I?  l  xn-MUDTSI  |H.  |N-1>  >  2002.2002.60 
2002  *<i*JDT$< |H. IN-1) 

13H.10 

*l«N«l, 

U3  Tu  500 
510  *<1«*X 
►1«F 

*<«*4*DEL4*SIGN 

IDxr-mDOTSIIN.  IN-1  >  >  2002.2002.60 
6 3  l“l*K-XAKAK>  100.100,520 
508  iri  |t»4AC-10>  709. 710.710 
71 J  ► 2,rl 
*<2**41 

702  ► l*r 

*<«  AdSUXK2*Fl-X8*F2)/<Fl-F2>  > 

|ri|KTMd-M530.529.53Q 
52V  H3UT»< IH, |N>  *Xk 

P  =  m  |UD'.  |K>  IN)*  (2.  •3,14161/ (X4*C> 

I DISH*10 

U*P«2 

l.*l 

^3  Tu  8000 

901  |F<|mTMB  -1)531.528. 531 
526  C*C-HTRBC 

l“tC-CH|N>520.520.546 
546  l“UN-3>532.533.533 

532  *<*  HOOTS! |N, |N> 
l»Ud  *10 

I3R*1 

KUK-XKNM)  500.500.61 

533  *<1*  HOOTS! |N. |N-2) 

*<  3«hQD  TS( |  M, |N) 

*<2*H00TS< IH, IN-1 ) 

Cl*  CM AX  1 -C 

Ol*  - ( XK3 -*4 1 ) / <  2 . *DELC  > 

U2  *  I  *  43  -2.»XH2  ♦XK1)/DELC‘>*2 

*<•  *41-  Dl*Cl  ♦<D2»Cl«(Cl-DELC))/2, 

1 ? I Rd*10 
l  3P*1 

l'<  XA-XKMAX )  500.500.61 
531  C*C-OELC 
l»Ud*l 

I* IC-CM|N)520.521.521 

521  IDIN-3)  522.523.523 

522  *<*4UDTS< IM, |N) ♦, 5* <R0OTS( IH. JN)-RDjTS! IH, IN-1 > ) 

1 N* 1**1 

1 3^*1 

in  xr-xkhax  >500.500 .520 

523  *<1«hOOTS! |M, IN-2) 

* <2 *  ROOTS!  |M, |N*1) 

*43*  RODTSUH.IN) 

C1*C"AX1-C 

CiAXl  aCHAXl-DElC 

Ul*  -!X43-X41c/<2,  •  UELC ) 

U2*  UK3  -2,*XH2  ♦  x Ki ) /  DELC  **2 
1 M* I N •  1 

*<*  X41  -01 *Cl  *(D2*Cl*<Cl-DEtC) )/2. 

13H«1 

|r (xr. xnMAX >500, 500.520 
520  ID IH-NOOE INNODE) >524,525.525 

524  I  MUD  *  l*'OD  *1 

n^mmuD  *  1 


1 1 1-6 


u«nuoe  (ihuo) 

‘XiRUOlSi  IH.il 

»<"|Ni,oliXO 

»<MA«il0.iROOTSt|H,l) 

1114ml 

CNAX1  icmax 

CaCNAX 

IDlSHilO 

lL*l 

1 JHPlS 

IN»1 

IXiNLAYER-1 
1 3! ALi f  THICx 
U3  TU  1000 
30 26  CiL'OAX-DElC 
IN*2 

lliNLAYER-1 

UHil 

U3  Tu  2000 

530  «i* ROOTS! IN, |N)i|CiPTRBC) 

■2>lmC 

U<  ■  HOOTS! IN. |N)  -XX 
Hi-tn2.Mll/0K 
Ci  C-PTRUC 

I3HH12.<3.141A/! I M2-M1 1/2. «M1) 
l  =  IH«H 

X<NAXiHDDTS(  IN,  IN)*10,KHDDTS(  I  *4 .  IN) -ROOTS!  t  .  IN-1  >  > 

x<hin..oi*roots;ih.in) 

hr  IN'  1020. C> ROOTS! I H , 1 N ) ,PEH I Ufi! I H , I N  )  .  R ,  T QRR 
1020  F3MMAT(l2l,Fl0.4.3X,Fl0.9,4X.F10  .4. 6X.no.  4. 17X.no,  4) 

U3  TU  S31 
61  UiC-HTHOC 

HR1NI  1031. C. ROOTS! IH, IN), PERIOD! IH, IN) 

03  TU  020 

025  HJNCH  2701. IRUN 
HiUSt  7 
HjAO  300,1  PUN 
03  056  i ■  i ,  i pun 

HXAD  2701 .MODEL  elH.L. TOTAL, C.PcRIOoIl.l) 

"s»0  1032,IU0ISRIJ),MDISPI J), J«1.L> 

HRINI  097. MODEL. 1H.C. PER  I  DO  1 1.1 1 .TOf Al 
os)  ) 3HHAT i ton-model  rti  ia  /,  oh  node  M/  ,ioh  phase  yeldcity  ro,3  /, 

1  7rt  hjr i bo  F«,J/  , 22H  SECTION  > H 1 C* Nf SSI XM >  p».3  ///  ) 

HRINI  2703 

2706  ) DHHAT  (  02H  DISPLACEMENTS  a*  LATEH  INTERFACES  NORMALIZED  TO  The/ 
151H  VERTICAL  DISPLACEMENT  at  The  SURFACE  I  INTERFACE  O)  /  ) 

HRINI  2704 

2700  iohoatiioh  interface  .ox.ish  hjhu.  disp.,ox.i2h  «eri.  dish.) 

U3  556  J10.L-1 

HRIN'  006  .J.UDISP! Jill.MDISPI J-l) 

006  )  3h 0 A T  (0X.13.11X.E13.6.9X.E13.6) 

HRINI  1024 

1024  I3HHAT  I19H  END  OF  COMPUTATION  ) 

H»USE 

enu 


INPUT  DATA  ANT  FORMAT  FOR  STRESS 


CARD 

FORMA 

DATA  DESCRIPTION 

1 

15 

No.  of  layers  in  model 

15 

No,  of  (e,k)  pairs  for  which  modal 

shape  is  to  be  computed 

2 

8F10.4 

Layer  parameters  stored  in  sequence 

Si  8 i  p  ,  d  ;  two  layers  to  a  card 

3 

E13.6 

c 

E13.6 

k 

E13.6 

Horizontal  to  vertical  surface  dis- 

u 

placement  ratio, 

w 

0 

u 

Repeat  card 

3  for  each  set 

of  c,k,  ~  . 

w 

o 


* 
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l*vrl :  l.inLhft 


IBB  05/160  BASIC  FORTEAB  tY  (E)  COBFILATIOB 


DATE:  66.  114 


R. OiM  1 

s.ooo  ’ 

5.00*  * 

5.00<'4 

S. PPOS 
S. OOP ft 
g.0007 

s.oom 
‘.OOP a 


R. OPW 
55 .  1011 

g.  m : 

•;.0i"  i 
5 .  .'p  14 
5,00  IS 

-a.0016 

5,0011 

5.0018 

5 , 0  0  1  *> 

5.0020 
5. 0021 

S. PP7? 
5,00.’  1 
<.?{'?  a 

5. CO'S 
5,  ''PIS 

5.  joi** 

^ ,  M"** 
<  m  v' » i 
S .  "  1 ' 
5.  >0  1 1 
51,  V  17 
5. 00  1  1 
'  .  M  1 4 

: .  ms 
■? .  2  *'  * 

s, Oft  t* 
5.  I'M* 

5.0'  10 

\  -'"JiG 
-*.  ‘  U  • 
5.  ’  -i  ■* 
*.  O'*  •*  ' 
5 . 0  1 4  4 
5.  OP4S 

s.p-m 


r  PRPGBAB  GTHESS.  .1.  tf .  DUMB  ntl-161 

DTHEBStO*  AtFAtSP<,BETARSO<,DR50)  .  > HORSOS,  LA B BDA«S0<,  AB«4 ,  4<, 

1  BO (SO) 

MEAL  LABSDA.BU 

DOn*LE  OBECISIOB  Al.BATtO.5tOa  TAOU , HBDO, BBHO.BEBSXG, BEBTAU, BBUB 1 , 
ICOSrS.COGQB.SIKPB.SIlQB.UDOT.UDOTB.HDOT.HDOTB.X.Y.PB.QB.BALPAB 
OOUM.E  PRECISIOB  DEI PPB . B ECIP P , DEIPQB , BECIFQ, BB BTiB,OB ’ , BBC2, AlPA 
DOUBLE  PBECISIOM  BBO, BOB. FEBtOD 

00ft B,  E  FBECtSIOB  C,  ALFBA.  BETA , B ATIA, K , 0.  OH, 0AEBAB,T*3TEB#C3QABE 
DODSLE  PBFCI510B  BEMfl ATl4< ,  OLDM?**< 

1‘»ft6  COBT IB*1F 

READ  (1 , M  LAY  EB5 , BCA3  55 

r  BCASFS  IS  THE  BDBBEB  OF  CASES  OF  C,K,BATIO  FOB  A  GtVEB  SET  OF 
C  LAYFB  F ABA B ETEBS .  ALFA.  BETA,  BBO,  D. 

S61  BEADS) . 1< (ALFARI<,BETlRt<, BHO  (t)  .D  (I)  ,  Y  *  1f  LAYEBS) 

BBITE(I.'i)  (ALFA  ( t)  .BETA  (I)  .1110(1)  .0(1)  ,  I «  1,  IA  YEB3) 

1  'OBBAT (HF10. 4) 

1  FOOBAT ( 2T5) 

EBT*R  *  LAYERS  -  1 
SO  200  L0-1.BCASES 
B  EADR 1 , 7<C,  K ,  BATtO 
C50ABE  *  C»C 

PERIOD*  h.20nRSAO72/lK*C< 

UDOTMATYO 
BOOT# 1.0 

OT.DBAT  S1<  •  UPO* 

OtDBAT  «2<  •  BDOT 

2  rORBAT  I  IE  1 1 . 6< 

HOTTER 1,fl<C, PFHIOD 

a  FaE3A?«1”1,SX, 1SBP»ASR-¥ELOCITY*,E20. 1 2,/, Si, 7BPEBI0D* . RX, E20. 12, 
1//, 10I.BHIBTFBFACE, 11,101, 1 2HH0B IZ . Dt SP  BI.IOI.IIBYEBT.DISP. ,/< 
B*0 

YRTT*  S3, 10<B, DDOT, BOOT 

10  ^OBBATRIOI,  2X,XS,**I,1X,6X,E20«12,6I,E20.  12< 

DO  100  B« 1 .LAYERS 
B,!*B<*BET.t  RR<  **2  •  H HORB< 

LABBDABB<*  RHO*R<  •  RAt*ARB<**2-2.0*BETARR<**2< 

100  COBTIBDE 
1701  5100*0. 

1202  TAt!U*0. 

1  -»ni  Qf.DB AT  •  SIGH 
1204  r,LDBA-  R4<  •  TAOU 
'MS  DO  20p  B  -  1,  KHITR 
1 70S  B 1  *  B  ♦  1 
17*»7  ALPHA  •  AL*AIB< 

17P«  ROB  •  SflOR«< 

1  ?Cn  3 AT! A  •  BETA«R< 

1710  DB  •  PRR< 

1711  X*DA*R ( (C/ALPHA)  **2  -  1 . ) 

1712  Y  *  DABS  ( (C/MTEA1  ••2-1.) 

1?n  ■»  ALr  AS  =>  D30BT (I) 

1214  BRETAB  -  D50BT  (Y) 

1?1S  GABBAB  •  2. *RB ATEA/C<* • 2 


mm  1216  P-  •  K  •  ^ ALFA B  •  OS 

mm  12W  3B  •  *  •  *»BSTA  B  •  DB 

mm  171*  tfiT’1  •  C  -  A  l p  H  A 

mm  1214  tf  (C* ALPHA)  22S.220.220 

’.70S 2  77"  C03"B  *  DCCS(PB) 

miSI  mfs  •  D5IB1PB< 

m"S*  AfASYH*  1 .0 

-.MS*  GO  TO  7J0 

-’.  VS  27S  DEtPPH  »  DEXPRPB< 

moe*T  RECIPP  •  I.O/DSXPPS 

7.205R  CCSOR  •  «r»EXPPR  *  RICI?P</2.0 

5.1PS1  SIRPB  •  EDEXPPH-RECIf>P</2.0 

7..  ■*'*«. 7  AFA5YB*-  1  •  0 

->10  CORTIBM* 

5.0''*?  TESTER  «  C  -  SATEA 

’.'*^*•1  IF  (C-BATPA)  2?3,717#217 

:.On‘.4  271  DEXPO«*DElFiOB< 

-.‘MS  PECIPO*  I.o/DEIPOB 

CC5QB  •  BDfXPOR  s  *!*C:P0</2.0 
ST RQB*  RD!IPOB-RECI®0</2. 0 
^.0 OS'  PTAB I H 1 . 0 

5 . 0 0 S  7  GO  TO  2^7 

'T  SIBOB  •  DSIRR2B< 

-."'/’I  C050B  =  OCOS(OS' 

.  ”  RTA5IB*1.0 

717  COB^IRnp 

TJ  G-1  •  GABBAB-1  . 

~.0D?C  RBC2  ’  FHO  (B)  •CSQAPE 

G.°G7f,  mi  A 5R  1 , 1  <  GABBAB*  C  )SPB  -GH  1  *C05QB 

-.007-»  111?  AHR1,?<*<r,si  •  5IRPB/BALFAB  E  GABBAB*  BBETAB  •  SIMQB* BTASYIK 

mi  mi,?'  ■  -rcospb  -  cosoh<  /  shc2 
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S.OG  »*'  BOD  7  BP  ■  BODRI  BP,  2< 
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“i.iini  do  ansi  kit  •  1,4 
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STORAGE  RAP 
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REL  AOR 
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REL  AOR 

RARR  TAG 
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000103 

PR 
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000110 
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r p* 
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RECIPP 
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0  0  M  4  o 
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000859 

APASYR 
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"  M  •:  v  «{ 

owisc 

rootrp 

000B60 

RCDRXY 

000869 

EXTERHAL 
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RARE 

REL  AOR 

RARE 

HU  AOR 

RARE 
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-  ••* 

POO*** 
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000883 
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00C99C 

«  *  ?  '  '  1 

099  VO 

oooooooo 

0008C4 

416937 Ed 51 121 B7p 
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INPUT  DATA  AND  FORMAT  FOR  INTEGRAL 


CARD 

FORMAT 

DATA  DESCRIPTION 

1 

110 

No.  of  layers  in  model 

F10.4 

Factor  for  singular  matrix  criteria 

in  curve  fitting  subroutine, 

generally  *•  .99. 

110 

degree +1  of  the  polynomial  fit  to 

“he  displacement  data 

T10 

Read  option;  1  if  read  every  set  of 

displacements;  0  if  read  every  other 

set 

2 

8F10.6 

Layer  parameters  stored  in  sequence 

o ,  B  ,  p  ,  d;  two  layers  to  a  card 

3 

20A4 

Text  identification  S  80  characters 

4 

14 

Model  identification  no. 

14 

Mode  no. 

14 

No.  of  layers 

F9.3 

Total  section  thickness 

F9.3 

Phase  velocity  in  km/s 

F9.3 

Period  in  sec 

5 

6E13.6 

Displacements  stored  in  pair  sequence 

vl»  w^,  U2»  . . .  ;3  pairs  to  a  card 

For  computations  involving  more  than  one  set  of  displacements  corresponding 
to  a  (c,T)  pair,  repeat  cards  4  i  5  for  each  set. 
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S  .  n  7  "i 
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Br AO  (7 

r  *i  7 ,  n  * 

•5,  '  n 

Pf  *■'  (1 

^ 70 

T»(\er 

■7  ,  *  "  1  0 

.*71 

no  "o 

<•.  '47 

:  IV 

T 

3  •  *7  *  4  1 

-  -  -r. 

3 . 0  r  -4  1 

**r  op 

S.  0  014*4 

Q*1P 

•J"Tc  • 

S.orjuE 

1 

P0"P  A“ 
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"0  1  "4 
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■»  ’ 
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: .  •'  o  *  j 

a 
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S . C  u  c 

HZ  *  0 
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S  .  0  o  e  1 
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*t.  0i>3  » 

•  •  «  7- 
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4  41  0 
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S.  0'*«,r 

r"*  A ,  * 

SA®»  • 

" 

0pT  ' 

''.O''*’’ 

H7  - 

r  2  *  n 

C. AO-  ^ 

or  i  ’n 

<: .  v  , » 

■*'  •  •* 

s. 

S  ■  n  7  «  1 

no  7  70 

a 

O  A  P  /.  > 

*'.77 

0*  t  * 

h-:t'  ( 

*4 

•  ■»  *  - 

S.n’*'  s 

7  70 

3.00*,*, 

2  7 1  *  •  , 

•.00^0 

IO.a 

'"OH-’H 

•  • 

o 

0*0  • 

vr-  i  «r  «  « 

-  a  r  1 

C  A  ’  !  - 

",  7 

t  1  •  ’ 

r .  0  :  V 

*>0  H»1 

0.007 > 
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0 , 0  "  » '* 

*»»  «f  /•  V 

0  .  0  7  0 

r  T  «  r» 

"AN®  •  xi.  rj-»  t  **»*;  t¥  THE  RGDEL 
P»N'"M  •  MO.  0®  HfACES  I*  RODEL,  t  RA*Y  R  1 
LAY’*-?:  t  mo  nr  **/•>!*  ACROSS  NRICR  U  AND  N  ARE  APERGX I  RATE  D 
5  r  VII  ■»•:■»  ,>r  jt,  ■»  *•»••**;  «s  ON  Mini  ’J  AND  H  APE  APPRQI l RATED  •  tAYrPD  F- 
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m  t. a y  r ^ 
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“1** ,  3’tVNT,  X,“ST-x,  ZZ,OPSTRNfOPRIR,  A,  BATE  A,  DC  DO, 

IS  fxn*#  7  A  Vf  ,  os:  V  00,1  A13D\,DiDD,2,’lllPPIG,  •),»!,  ALPHA,  DCOALP,  At®  A, 

ttGTinPI.,3  7.  ?*»S#*A-:?OR 
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INPUT  DATA  AND  FORMAT  FOR  TRAVEL 


Card  Format 

1  14 

2  20A4 

3  14 
14 

F8.4 

F8.4 

F8.4 


F8.4 

4  10F8.4 

5  10F8.4 


to  be  computed 

Identification  for  models  £  80  characters 


NLAYER  ■  no.  o.  velocity  layers  in  model 
NR;  depths  of  penetration  for  which  T-D 
values  are  to  be  computed  *  (NR)*(DELR) 
DELR  *  Interval  in  (km)  for  numerical 


integration 

DEPTH  ■  maximum  depth  (km)  of  penetration 

£  depth  of  last  velocity  value 

EPSLON  «  distance  (km)  from  maximum  depth 

to  begin  using  modified  trapezoidal  rule 

of  integration;  usually  .3km 

DELTA  ■  distance  from  maximum  depth  at 

which  integration  is  stopped;  usually  .01km 

Velocities  in  km/sec 


Depths  (km)  from  surface  of  given  velocities 
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0 I  HENS  |  UN  T(200).UEi<200)>B(2QO).VEl<200). 2(900). ANGIE(200),V(900> 
1. 02(000). T£xT(20) 

O3J0LE  PRECISION  T, DEL.R.YEL. 2. v.0i.DElR.0EPTH.5PS.0ELTA.ni.R2.R3. 

7  0;  I  .  ».8.C.D  =  Ci  Hl.RO.XR.P.  P2.SU8.SUH1.SUH2.  ETA.  fact.  Mc’l  .F»CT2. 
2IE»T.YVl,00.D..ll.DEPTN?.fACT3,FACT(i.FACT3.fACT4,TA0.H0P,SAD.Al.A2 
3*il. 82. Cl .C2 

U3UBLE  PRECISION  2MaT(4,4),HDUHNY,YVEC(4),AVEC' «),0 
•10  E3aNAT(4E20.16) 

M-A0  (1.100)  NMUN 

03  300  IRON  •  1,  NOUN 

MEAD  (1.400)  I’EXTd  ).IH.20> 

"EAC  (1.100)  NLATER.NO.OELN.OEOlH.EOStON.OEtTA 
MjAO  (  1.200)  <  V  EC  < 1 ) * lll.NLAYEM) 

MjAO  (1.200)  (N< ) ). 1>1,NLATEN) 

100  F3MMAT  (214.4F0.4) 

200  F)H1A!  (10F6.4) 

4 1 1 )  sq 
I«1 

1*  1.1.1 

4( I  )«2( l*l>*DELM 
I*  U.l)-OEPtH)  19.19,20 
20  Oil 
All 
111 

IJ"P*1 
7  Jlj.j 

NRMiJ 

42  L3NTINUE 

03  no  L  1  •  1 . 4 
N0U4HT  l  NR4  -  4  •  LI 
DJM4»  ■  R1NDUMMY ) 

«3ITt  (3.0000)  DUMMY 
2m*T<l1.1>  «  1. 

03  700  L J* 2. 4 

44*I(lI.LJ>  •  UOMMYi.(LJ-t) 

7Q7  CONTINUE 
710  CONTINUE 

NOJMMY  ■  4 
•  0.001 

CALL  SJR4I2MAT.NDUMMY.ERS) 

03  710  I.Jil.4 

NQOMMY  t  NRH  -  4  •  LJ 

YYtytlJ)  •  YEL(NDOMMY) 

713  L'3-iTIMUE 

03  720  L 1 1 1 . 4 
UJH4T  I  o, 

03  720  LJil.4 

0J«4»  •  UJMNY  •  2m»T(LI.LJ)iO*EC(LJ)  , 

720  L'3«IhUt 

A  YE  C  <  L I )  «  DUMMY 
723  LOnTinuE 

A  ■  AYEC ( t ) 

8  •  »YEC(2> 

C  •  AOEC(J) 

0  •  AYEC ( 4  > 

(.3  Tu(31.40).IJMP 

31  V(l)  i  A  ♦  0.2(1)  ♦  C«2( I >**2  *  0*2 ( 1 ) i#3 
lll.l 

lt(2H)-N(J))  31.31.30 
30  1-*(2(1)-UEPTM)33.33.34 

33  (-'(Ml  J.3)-0?PTH)  7.7,32 

32  X(J)iOEPTm 
I.J  Tu  31 

34  0-PTmi*e371 ,-OEPTm 
1  ■  0 

37  1.1*1 

21  I  ME371  ,-2(  1  ) 

1  - <  2 ( I )-DEPTmi)36. 34.30 
36  1**1 

83*6371 .0 

A««MH 

hC<  OEPT*/(Xfi»OElft) 

*«1 

MPl  s  NW  ♦  1 

1 nC°1  >  INC  <  X 
M*  I  4 

U3  8  J  >  NHPl .  JKCPl .N& 

U-tl«OEL« 

1 J«a«2 

U2(Nj  i  6371. -2(J) 

M  t  XK»OEH 

I  7  I637t .0-/IJI-KINRR  >>  42.42,4* 

45  *M«  *  NNM  ♦  3 
U D  Tu  42 

>*2«p»»»2 

If  M 40/V( 1 M*»2-F2 >  4Q00.40*,4U!> 

407  LSnTlNUE 

Sj"*USQHT< HD/¥(1)>»»2-P2> 

\“ <  J-2>11.11.10 

11  >J«1*0. 

*JM2*0. 

3U  10  12 

iu  ij«l*l ./< *0«SUH> 

sj«2*SUH1*(R0/  V(1M«*2 
k*2 

24  feU**OU/XH) 

If  1 1 T* -P  )  4000.240.240 
240  CONTINUE 

►  *CT  «ubORTUET*-P)*(feTA«PM 

►  actiu./ ifact*2iri  ) 

UCT2*  FACU  *(2U)/V  (K))«2 

5j«i*bUMl*2.*FACTl 
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6jrt2*SUM2*2.  *FACT2 

*  *4*1 

|f(ZU)-2t  J>-EPSL0N)22.23*24 
2 2  ij:>T*Z(J)*EDSLON-,01 
IM2«0-TESI )12. 23,23 
2.1  U^Thi*^  (  J)  *EP5lCN 
U=PT«i  *6371 .-DtPTMl 

y^l«A*DfcPTHl*(B*DfcPTHl»<C*DkPTHl*OM 
USPTHl  *6371 , -DEPTHl 
U0*c^SLON 
tTA*UEPlHl/VVl 
1  •*  tfcTA-P)  <000.230.230 
230  CONTINUE 

►  ACT*DSt)»T((6T*-P)«(ETA*P»  ) 

► ACTl*i./(FACT*DEPTHl) 

FaCW*  FACT1«(DEPTH1/VV1  >**2 
6J*1  *S>UM1«UELN  ♦  FACT1«DEL1 
6J*2  *  SUN2*DELR  ♦  Fact  2*DEL1 
OD  Tu  21 

12  USL1*-2(J)  -  EPSlON  *Z(K-1» 

03  To  23 
21  1*1 
R3**0. 

03  *tPSLDN 
U;HTH2«2( J)*D0 

14  H3**KD*M, 

1*(  1-1)16.16.17 
16  Uir'THl*DEPTH2 

Uir'THi  *63  7  1 ,  -DEPTH  1 
Vyl*A*t  PTNi#(H.DfcPTHi#(C*D6PTMi*D)  l 
U=PTmi  *6371 ,-DEPTHl 
fcl AxuEPTHl/ VV  1 
1-  (tTA-P)  4000.160.160 
160  C  J**  T  1  NUfc 

►  ACf*DSUR!(  <ETA-P)*(ETA*PM 

►  ACT1*1./(FACT*0EPTH1) 

►  ACT 2*F act l*(DtPTMl/VVl)#*2 
1*10 

03  T u  l» 

1/  FaCT6  *  facti 
HCT6  *  FACT2 
► ACTi*FACT3 

F  AC T 2 *F AC T4 

id  U;»Mh2*2<j)  *00/(2.  ••PD*) 

US* Tm2  *6371 . -DEPTH2 
Y''2«A*D£PTm2«  (8*DEPTH2*(C*DtPTM2«DM 
U •  r' T « 2  *6371  .-utPTM2 
fc|A*UtPTM2/¥V2 
I?  (tTA-P)  4000. 160. ISO 
160  CONTINUE 

FaCT*DSORT( (ETA-P)*(ETA*P>  ) 

► ACT3*1./(FACT*DEPTH2) 

F aCT4  «F  AC  T  3* ( DEPTh2/VV2 )  *#2 
6JNl*SUNl»DD*2.*(FACTl*FACT3)/(2.»»PON) 

6  j«2*SUM2*Di)*2  .  •(FACT2*FACT4)/<2.**j0R) 

1*1DU/(24«*s»0o)-.03  H5.15.14 

15  ir<OU/2.**POM-  DEL T A ) 60 • 60 . 6 1 
61  IAU  *  00/(2. **POW) 

».0P  «  (V( J)-0V2)/TAD 
62  MCM  sFACTl 
fACT6  *  FACT2 

►  ACTl  «  F  AC  T  3 
-AC‘2  «  fACT4 
TAO  *  TAD/2. 

0  =  HTpii  s  2(  J)*TA0 

640  *  6LDP*TaD 
*Vi  *  V ( J) -SKD 

►ACT  «  2»*P*(V(J)*TAD  ♦  2<J>*b*0>  /( V( J)**2-2 ,*V( J)*SKD ) 

►ACT3  «  1./(DSORT(FaCT)*OEPTH1) 

•ACM  s  r  ACT3*  (  OEPThI/V  VI  )**2 
6J«1  *  SUM l  *TAD#<FACT1  •FaCT3>»2. 

6J«2  «  SUM2  ♦  TaD*(FaCT2  *FACM)  *2. 

!► ( I AD-DELTA)60 .60.62 
60  Ai  j\  AC  T  3 
A2  t  TACT4 

61  *  facti 
32  *F AC  T  2 
;i*) acts 
C2*F ACT  6 

ePND**  TaO#(FaCT4*3«A2-5.#B2/2.*.5»C2)/SUM2 

*3 1  Tt ( 3 • 1050  )  ERROR 

6JH1*  Ta0«(FAcT3a3.#A1-5.«61/2.a.5*C1)*SUR1 

*J*2«  TAD*(FACT4*3.*A2-5.*B2/2.*.5*C2) *SU*2 

U=k(N)*P«SUMl#6371. 

\  ( n ) *SUM2 

6  I  fli  I U  *  P*V(1)/RD 

C36IU*  SORT ( 1 ,  .  S I N |0**2 ) 

ANuLb(N)  *  ATAN(SINlO/CDSI0)*l60./i.l416 

NrN*i 

U3  TU  6 
400J  CONTINUE 

U;L(N)  0. 

HR)  *  0. 

A4<sLb(N>  *  0. 

*  *  L  ♦  1 
6  C  jt> T  I  SUfc 

N«N-1 

“41 T  t ( 3 . 40 1 )  1 1  EX  T ( I  ) ,  1*1,20) 

»R i Tb ( 3 • 1000 ) 

1000  FORMAT  (  6X.10M  DEPTH!**)  .  10* . 1§hTElDC I TYCJCM/SL,.)  ) 

M^|Tt (3.1050  )  (2( J) ,V( J),J«1, l I ) 

105 0  FORMAT  (2X.E13.6.6X.E13.6) 

m4 1  lb ( 3 • 401 )  ( IE  XT (  I)  •  1*1.20) 

M-UTt  (3.2000) 

2000  F3MMAT  ( 2 X c 1 6h  TRAVEL  TI*E  (SEO  .54.11H  OELTa  <KH>,5X,21m  OEPTh  0 
IF  HcNETHATION  ,5X  ,1VM  ANGLE  uf  incidence  ) 

M4lTt(3, 2050 )  (T( j) ,OEk< J) .  D^ ( J  > .ANiLE(J) .J*i .N) 

2050  F jnnaT  (  5X.F6.2.13X,F6.1.14X»F7,1.16X.F7.2) 

401  F3H4AT  ( 1 Hi  //// , 20A4 , ///  ) 

400  F  3HNAT  ( 20A4  ) 

300  CONTINUE 
61  JP 
b  mU 
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INPUT  DATA  AND  FORMAT  FCR  VARGRAV 


CARD 

FORMAT 

DATA  DESCRIPTION 

1 

13 

No.  stations 

F10.5 

Station  spacing  in  km,  Staspa 

F6.1 

Depth  to  mass  sheet  in  km  ,  Depth 

13 

No.  of  Density  layers  S,  50 

13 

No.  of  Density  profiles  £  100 

F10.5 

Minimum  $  value 

F10.5 

Delta  z  spacing  for  integration 

F10.5 

Subcrustal  Density 

^  min  ~ 


[exp(  Depth/Staspa)  ~  1.]  *  Depth/Staspa 

2 


No.  of  t's  to  be  used 
2 


2 

3 

4 


1QF8.2  Gravity  values  in  milligals 

4012  No.  of  Density  profile  to  be  used  at 

each  station 

16F5.2  Density  values  stored  by  profile 


III-19 


UjHgNSfON  GMaV'  100  >»DENSTY (50.50 ) .*C0RR(  t 00) ,0EP0 IF <100 > .PH | (50 >. ! 
iH*UFt 100  ) 

WrAO  100. *OSTA.STaSPA, DEPTH,  IlifcNS.MKOF  ,PH|*In,0El2  .subcru 

100  ► 3"4 AT ( 13»Fl0.5.Fto.l»2!3»3FlOr5) 

«  E  A  L  200*  (G**v<  I  > .  |*i,NOi>TA) 

200  »■  OrtiAT( 10T8.2) 

H=AO  150, ( JPHOF ( ! ),!«!. NObTA) 

150  )  3*4*1(4012) 

MS*0  300,  * (OtNSTV( i,JI, J»i, !Ofc*S>. !«l,NPROF) 
iou  *  0H1KT< 16F5.2) 

UD  1492  1*1, NPNOF 

U3  1492  J«l,  | OEMS 

Or N5 l Y  ( 1 , J) «  SU8CRU-0ENS7Y  (t.J) 

1492  CjMPiUfe 

t*HOUO  =  i -feXP<  3, 141**  OEPTH/S?AS^A)-i.  )«DEPTh/STaSPA 
fc*PEv'4«  *EAP(4, 1416*  DEPTh/STA1>PA)-1,)*0EPTh/STA  SP* 

A  1  *  0  - 

I  *  1 

^ AM T  *<0EpTh**2)/(STaSPA**2> 
a  l‘t  (I  ♦  1 )  /2  -  1/2)1, 1.2 

1  P4l( I )tfeXP000/(PART*A!*»2) 

P=l!NI  10l,P«I(|) 

101  )3M4*T  ( F 10.5) 

U3  Tu  3 

2  P4l( !>«fcXP6VR/(PART  ♦  # I »»2  > 

102  F  DM4»nFlC,5) 

3  I* < AhS(PHI  (  I ) ) -PM! M|N >4,4,5 

9  1*1*1 

*  i  *  A  J  ♦  1 , 0 
G3  Tu  6 

4  111  *  2*1-1 
Ik  «  1*1 

1X2*! 

J  *  I  X 

mid  1*2. 1X2 
^11  J)sPh|<  I  ) 

J*  J.  1 

id  13NT1VUE 

PH1UX2  >«PM|<1> 

Jxl  *1 

I *2*1*2  * 

03  ;  !?1,|X2 
r* l» !)*PM!( J) 

J«  J-i 

1  U3<*T  i  9Uk 

400.0EPTm.STaSPa 

400  ►  3*4aT(2)M-CjHPUTED  Ph*  FAC  TOtJS/lOM  OEPTKf  K* )  ,F6 . 1 , 2  A ,  20M  STATION 
i*a*ClMG(*N>.Fi0.5> 

pm*  <oi.(pHi(n,i«i,i*i) 

401  ►  3*4AT(2X.E13.4> 

1*1-1  * 

1  **<  *N *  I  * -1 

1  X  2  *  A  - 1 
D3*>  l«i.!X2 

*:jrh( I >  *  o . 

03  9  J«i,|x 

*:U«H(|)«  GRAV(  J)*PHI(K)  «XCO«iN<|) 

A  *  A  ♦  1 

V  COM  INUfc 

500  ► 3"4*T(2X ,F1P.*> 

I  X  *  I  *  *1 
A*  N  «N-1 

0  continue 

lx3«N05T  *-  |M«i 
Jti X2«l 
UJ  iu  |«1, |X3 

*:jmk( j)*o. 

A*  1 

1)3  11  L  *  1 ,  1*1 

*-URH( j)*GR*V(X )*PM| CL)  ♦XCORR( J ) 

*  *  A  ♦  1 

U  CDwTlNUE 
600  > 3HMKT(2fc.FiO,2) 

J*  J*i 

10  C0M1*U£ 
lxJ«lX3*l 
n« p»0b  T  A-  I  *2 
03  12  !«!*3.N 

*  2  I 

»CjH*( J>*3 . 

03  10  L  *  1 »  H 

*:oMA(  J)«QH*V( A )*PM|(L >  *XC0RP( J) 

A  *  A  ♦  1 

13  CONTINUE 

?C0  ► 3«4AT(2X.Fi0.2) 

J*  J  *  1 
MzA-l 

12  C3NT1NUE 

f  ACT *2. *(3.14 10* *2) *6. 6 7 
L‘J  14  1«1,N0STa 
*;uM*(|)«  XCOhH( 1 i/F ACT 

14  continue 

03  JO  JU.NOSTa 

5  JrfMJ  . 

AiA  V*  XCOKR ( J ) 

»’( AHAX-O. >50,51.52 

50 

03  Tu  53 

52  V;r**i.O 
U3  TU  53 

51  U?PO »  F ( J ) *0f PTM 
1*3  Tu  30 

53  !«1P"0F(J) 

CrHOlFCJ)*  0. 

A  *  0 
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40  K»A*1 

U-NS  ■  DtNSTYJKK) 

U=NS:*DEN5Tt<  1  .  X -1 ) 

K  <K-IDEN$>43.43.44 
*  D?**D1F (  J)  *QERDI F  ( J )  *DELl 

agH*  SUN*DENS*UEL4 
M  AbS  I  SUN)- A85UPAX  >)  40,41,42 
42  |r(*-l>»4. 55-54 
55  U;NS1*DENS 

54  P.AgSUNAXl-SUH^DCNS^DEH 

US^d I f I J> iDtPT 4-SGh  ' I DE*Dlf  U> ♦»'Dfc YSi-UELZ > 

GO  TU  30 

41  UjFDin  JMDEPTM-SGN*  DEPDIF(J) 

U3  To  30 

44  y.*$  .  DENSTYU.K-l) 

Ur«Sl  =*DfcNST  Y  (  l.K-l) 

IWrt-1 
CD  TO  43 
30  CONTINUE 
♦Ml*)  2000 

2000  U«NATtiH-,43X.27H  THICKNESS  OF  C»UST*L  LaY(R  /) 

P4HI  3000 

3000  ►3hNAT(17H  DENSITY  PROFILES  //  ?X.1QH  0EPTH( KH ) , 5x ,  9h  OENSlTY  > 
l‘;PTn*  -06LZ 
U3  35Q0  1*1.  NPROF 
DEPTH.  DEPTH  ♦  DEU 
p-UNI  360  0. I DENST Y  ( 1, J>, J»l, 1  DENS ) 

3500  CONTINUE 

3600  ) OHNAT  (  10F9.3) 

♦MINI  196b 

196b  f 3«NAT(bH-STAT 10N.5X.19H  G»UV 1 T Y ( N|WU GALS) . 2X , iQH  DEPThCKHI. 

1 2 x , 2 1 h  DENSITY  PROFILE  USED) 

CD  1/  1*1 , NOS  TA 

'MINI  4000. l.GHAVt  1  ),  DEPDIF1  1 ) . 1PHJF{ 1 ) 

4 uOO  *3*14? 1 2*. l3.13X.tM.2,9X,F6,2,13X, 13) 

17  CDNTlNUfc 
PAUSE  / 
tNU 
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INPUT  DATA  AND  FORMAT  FOR  PRESID 


FORMAT 

13 

16F5.1 

20F4.1 

13 

F6.1 

F6.1 

18F4.1 

18F4.1 

16F5.I 

12F6.1 


DATA  DESCRIPTION 

No.  of  stations  at  which  residuals  will 
be  computed 

Jef frevs-Bullen  P-travel  times  in  sec 

Latitude  corrections  for  converting  to 

geocentric  coordinates 

No.  of  epicenter  locations  for  each 

station 

Station  latitude  (geocentric) 

Station  longitude 

USGS  origin  times  for  earthquakes  stored 

in  sequence  as  HH:MM:SS.S 

Record  arrival  times  and  store  in  sequence 

as  HH:MM:SS.S 

USGS  focal  depths  in  km 

USGS  latitude  anJ  longitude  of  epicenters 

stored  in  pair  sequence  Lat^,  Long ^ ,  Lat2» 

Long 2 i  . . • 

Sign  Convention 
East  Long  - 
West  Long  + 

North  Lat  + 

South  Lat  - 
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u("E*<>10N  UOAL0C(2Q0) .ORT 1H(3Q0 )  .  Ue>*TH  <  ICO  >  .  $?  A  i  i  H  (  3QO  1  . RE  S !  0<  100 
1  * ,  M  >E*<  100  >  .  DELTA  tlOO  >  .a2»KuU1QQ  >,  9  <14 . 103  > , IOCOEl  <  1 00  > .  GEOC<  9 
11)  .i*edUA(20U) 

«s*0  100.  UHSTa 

lOU  E3*<4AT(13> 

RjAu  150, C [PC! .j). I ■1,14), J«l.lo3> 

150  E OhiaT ( 16^5. 1 ) 

h?aD  U|2,(UEUC< J).J>1.91) 

.11,?  E  3"4AT  (20^4.1) 
w-  s:  I  STATION  IKUwX 

U3  1  1«1.U«STa 

"tAQ  <00.  lihooa 

200  OH^aT(I3) 

U«?*Ll*0UA 

*  c  A  0  700  .STALAT.  STaLRG 
700  E  3*HAT<F0.1.F6.1) 

U4  *  3«  tMOUA 
H;A0  600.  <OHTIH( jl . J«l, 1X2) 
oo  fdhuat  ti#f4*i) 

«  =  AD  400.  <i>TAl  I H  (  J)  ,J»1,  1X2) 

<  u  0  » 3H4AT(iaF4.1 ) 

*iAQ  500. <0EPTh< J) ,  J* 1 , L I MQUA > 

500  E  3* 4AT ( 1 6F5 , 1 ) 

HiAO  300. (OUAIOC< J). J«l. IX) 

300  Y  3H1aT( 12F6.1 ) 

03  1111  J«l. IX. 2 
Ka2 

AWAT  *1,0 

502  )’«auAinc(j)-ALAT)565. 50 1,501 
5U 1  A. AT  a  A  .T  ♦l.O 
A  bA  *1 

GD  Tu  90? 

663  COMb;EQCII<>  -(G£0C(K)-6£0ClK-l  )  )»UL»t-UU»LOCt  Jl) 

1111  =  OJALOCU>  •  CORN 

U3  2  J»l,l*.2 
I  MJ»sj/2»l 

3Tx»n3  ■  S7»L»I*3. 1416/1,0. 

UJANU  »U£3C*(JI*1. 1416/1.0. 

•*  JL  ANG  •  3TALNG  -  aU.LOCCJBl) 
i-  (I-3L4N3-160.  13.3.4 
4  ''a  .»>lG«  ‘"O.ARG  -360. 

3  HDlA.bJ  iPJLANGb3. 1416/1,0. 

C330E.  B  Si  'll  STaANG  l.iim  UUaNG  I  .C06I  AANG>«CO!i(OUANG)«CO!><POlANG> 
^l>0tL<SO"Bl(l.-C0!>O£t**2l 

IKOXIb.TARISINDEL/COSOElI  •1.0. /3. 1416 
uHrA(|N0»>*  DfeL  T  A  t  [  NO  x  )  •  00  .  -S1GH( /U.  .OElTK  1H0«'  ) 

3 1  «5t  T bC03(  JLXRG  )*S|NIPQL»HGl  /»  1*<UE. 

CD.Jt  I  «<S|NIOU»RG>*CU6(ST*X»,GI-CU6l  JUAR5  I  »S  I  * «  S 1  X  A*r,  >  *C0S  I  POl‘NG  >  > 
1/51 MOcL 

AjIVb.bMSiRBED 

»• i txB4r,vi bsin/cosbei i»ieo./3. 1416 
B;lIXB  ti£TTX  .  00.-SIGNI9U..BEITA) 

I/IjIXBET  -0.162. .4 ,»3 
a*  iri;uit)£T  -o .  16. *63. BO 
63  B£ I T A  lias. 

93  A/IIUTI  1N0XU360.-BETTA 
uD  Tu  81 

62  Aj | AuTl 1N0« l«Bt TTA 
91  LJ.Ii'lUt 

2  CD  bI 1 9Ut 

-'ll 

JD  12  jxl.L'nbUA 
A  *2 

Ut**  s  33.0 

1 r ' Ut _  T  A ( j  )  - 1 02 . >  13,13.14 
13  i-  l  ufc  H  C J)*DEP)i5.15.16 
1  ft  U-h  *  OtP  •  63.0 

ft  rA  •  1 

03  ti>  13 

15  l*2 

U|>T  *  1,0 

1*  l-«0tjA(  J)-DI  ST  >17,17, 18 

le  0  I o  T  *  0  I bT  •  1.0 

lsL*1 
.3  TO  19 

1/  J 1 E  1  ■  PIK.D*  (P(R.L>  -P<A.w-l  >  )»tUlST-DEwTA(  j>  » 

UJE  2  «E(d-l,l  -|E(H.L)  -PCA-1 .1-1 )  )  •  C  D 1 3  T  -0ELTa<J)) 

UIDfcPTKt J)-  .0)  90.90.91 

90  '-tc3>4  •  0If*  -  (  0 1  F  1»L  I F2  )  •  il)EP-OrPTM(  J)  )/33.0 
«3  TJ  92 

91  ‘hc3H  «  Ul^l  -  T  0  I  T 1  -0IF2»»<0i-P-0fc‘>>H(  J)  J/63.Q 

*2  *T i 9t*0RT I«d^)»36C0 . •  UN »  I H  ( M  •  1 )  • 6Q . ~ORT I H ( h • 2 ) 

vTItfc  a  s  T  A  f  Ih  I  H  )  •  36C0  .  •  S  T  A  T  I  H<  NM  '  •  6U  .  »S\  AT  I H  (  H»2  » 

N|NO 

1-13' I  HE  -  P5562. 6 >20.4000.4030 
400-  I  " « i « I  HE  -  0J552.6)  21.20.20 
21  ifht  a  STl'dE  •  06400. 

2u  MrsIuJj>*  ST -OTJHE  -  THfcOtH 
03  TO  12 

C  BRANCH  OF  RESIDUAL  TJWE 

14  H:b|u(j)>0. 

"j'O 

12  CONTINUE 

c  ORjfH  EPICENTERS  b4  OITTaNCE 

*i  ieh(H«oelTa<;  ) 

*i  >e«;2)«oelTa( /) 

i:wOFim»i 

v.Jo0tL»2)  *2 

U3  60  ft*3.  LlMUUA 

P*1 

Lai 

50  l QclTaTX  )-XITem(U  >22.22.23 
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23  Ls«.*l 

IT  <L-K>24,25,25 

24  H*«*i 

‘33  T U  3J 

25  X|  IfeHlL  >«OELTa<K) 

L3CUfcL< H*l)  *K 

<•3  Tu  80 
22  IU*« 

*»i  4j  > eh<  ixn*  xi  lent  ixi-i ) 

LDCOfclUXlMLOCOEUUt-l) 

I xl«l xi-i 
l?t IX1-MJ27.27.51 
2?  M  l*H< !  U  >*061U<K  ) 

L DC Ot L <  I  XI) «K 
BO  CONTINUE 
*M.Nl  800 

8«0  *■  3«NaT(*5h-  4ECO»OING  SUTjON  CPQR"  /26*,9H  LA  T  I  TUqE  .  AX  ,  10H  LONG  I 

2 1 JUE  > 

R9|N«  900.  5TALAT,STALNG 
9U0  ►  D*4ATC27X.r f>  1 , 9  X  . Fft  .  1  > 

K9INI  1000 

1000  ►  3***T  < /8BM  gPlCENUH  PARAMETERS  ORDERED  0*  OISTaNCE  -aZIHUTh  HEaS 

UAfeO  COUNTER-CLQCAHl SB  FROH  NOR«H  //) 

^INi  2000 

2000  T3*NAT<3X,9h  DISTANCE,  JX.9R  L A  I  I T UU r . 3X , 10M  LONG  I T  UqE  ,  3X  ,  6«  OcPTh, 
1 O«*0w  AZIhUTm.3x.33m  P-NESlDUAL(DBicNV£n-TnEORETlCALM 
03  6U  K>1,1X,2 
UltK/2»l 

L  sL.lt OEL  (  *  xi  ) 

R*  2*L-1 

^hl  3000,  XI  TEN  (  IX1),0UAL0C<N>,UUA.0C««*1) ,  DEPTr  l  ),  aZiMUUl  ),HES 
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APPENDIX  IV 


CURVES  OF  PARTICLE  DISPLACEMENT 

CURVES  OF  PARTIAL  DERIVATIVES  OF  PHASE  VELOCITY  WITH  RESPECT 

TO  LAYER  PARAMETERS 
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-1.0  l.o 


HORIZONTAL  AND  VERTICAL  PARTICLE  AMPLITUDES 
NORMALIZED  TO  THE  VERTICAL  AMPLITUDE  AT  THE 
SURFACE  FOR  FUNDAMENTAL  RAYLEIGH  MODE, 
GUTENBERG  -  BIRCH  n  MODEL. 


Fig.  IV -6 


FIRST  HIGHER  RAYLEIGH  MODE  PARTICLE  AMPLITUDES 
NORMALIZED  TO  THE  VERTICLE  AMPLITUDE  AT  THE 
SURFACE  FOR  THE  GUTENBERG- BIRCH  II  MODEL. 


Fig.  IV-7 


FIRST  HIGHER  RAYLEIGH  MODE  PARTICLE  AMPLITUDES 
NORMALIZED  TO  THE  VERTICAL  AMPLITUDE  AT  THE 
SURFACE  FOR  THE  GUTENBERG  -  BIRCH  II  MODEL. 


Fig.  IV  *8 
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Fig.  IV- 9 


THIRD  HIGHER  RAYLEIGH  MODE  NORMALIZED  TO  THE 
VERTICAL  AMPLITUDE  AT  THE  SURFACE  FOR 
THE  GUTENBERG -BIRCH  II  MODEL. 


Fig.  IV*  10 


